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ABSTRACT 

We report on the discovery of a rapidly co-rotating stellar and gas component in 
the nucleus of the shell elliptical NGC 2865. The stellar component extends ~ O.blh^QQ 
kpc along the major axis, and shows depressed velocity dispersion and absorption line 
profiles skewed in the opposite sense to the mean velocity. Associated with it is a young 
stellar population with enhanced H/3, lowered Mg and same Fe indices relative to the 
underlying elliptical. Its recent star formation history is constrained by considering 
"bulge+burst" models under 4 physically motivated scenarios, using evolutionary pop- 
ulation synthesis. Scenarios in which the nuclear component is formed over a Hubble 
time or recently from continuous gas inflow are ruled out. 

A recent starburst can satisfy observational constraints only if its population has 
metallicity 2.5-6.3 times that of the bulge. The nuclear iron to magnesium index ratio 
can be explained by a temperature effect in the atmospheres of stars at main sequence 
turnoffs between A3 and F4, during which the Fe indices of the burst population are 
high enough to compensate for dilution effects. It is therefore possible to modify line- 
index ratios (and hence the inferred abundance ratio) simply by the presence of a young 
population with the same abundance. The high metallicity requirement suggests self- 
enrichment, and burst duration longer than SN II feedback timescale. No solution exists 
for bursts longer than 0.4 Gyr. Burst age estimates of 0.4-1.7 Gyr are larger than that 
for the shells (0.24 Gyr) assuming phase-wrapping. 

No starburst is required if the nuclear component is composed of stars with Fe 
abundance enhanced by ~ 0.08 dex relative to the underlying elliptical, which are 
accreted by an event which truncated the star formation. This relies on the abundance 
differences between giant ellipticals and spirals. The age-estimates of 0.1-0.4 Gyr in 
this scenario are in closer agreement with those for phase-wrapped shells. 

Our results argue for a gas-rich accretion or merger origin for the shells and kine- 
matic subcomponent in NGC 2865. Arguments based on stellar populations and gas 
dynamics suggest that one of the progenitors is likely a Sb or Sc spiral. We demonstrate 
that despite the age and metallicity degeneracy of the underlying elliptical, the age and 
metallicity of the kinematic subcomponent can be constrained. This work strengthens 
the link between KDCs and shells, and demonstrates that a KDC can be formed from 
a late merger. 

Key words: galaxies: interactions - galaxies: internal motions - galaxies: nuclei. - 
galaxies: individual. 



1 INTRODUCTION 

Observational and theoretical arguments suggest that shells 
in elliptical galaxies and kinematically distinct cores (KDC)s 



may have a common origin. Observationally, about 15-20% 
of shell galaxies have nuclear post starburst (E-l-A) spectra, 
suggesting that shell formation process may involve some 
accretion or merger of gas-rich material with associated nu- 
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clear star formation (Carter et al. 1988). In some ellipticals 



the onset of the kinematic peculiarity is associated wi th an 



enhanc ement in the metallic absorption-line str engths (Ben- 
der & Surma 1992, CaroUo & Danziger 1994), which may 
suggest that KDCs are products of early gas-rich mergers 
with an associated metallicity enhan cement. The s trongest 
association comes from the finding by Forbes (1992) that all 



of the 9 well established KDCs and a further 4 out of the 6 
"possible KDCs" possess shells. 

The models proposed for the formation of shells and 
KDCs are very si milar, and ma ny involve galaxy interac- 
tions or mergers, ^uinn (1984) considered mostly radial 



collisions of low-mass disks with a spherical potential, and 
found that shells could be formed through the process of 
phase wrapping. His ideas w ere further developed for low- 
mass spheroidal comp anions (Dupraz fc Combes 1986 ), and 
for non-radial orbits (Hernquist & Quinn 1988). Hernquist 



and Quinn further showed that mere mass-transfer dur- 
ing parabolic encounters are capable of producing shells, 
demonstrating that a merger is not necessary in the phase- 

wrapping model. 

Hernquist & Weil (1992) extended the models in which 
a dwarf is captured using a hybrid N-body/hydrodynamic 
code, tracking gas as well as stars from the captured galaxy. 
They concluded that, while the stars formed the shell struc- 
tures, gas would settle into a compact disc or ring in the 
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Figure 1. R-band image of NGC 2865, showing the bright, 
chaotic shells. 



When seen from different angles, the shells are either inter- 
leaved or all-round, providing a natural exp lanation for Type 
1 (alig ned) and Type 2 (all round) shells. [Hau fc Thomson 



centre i |)f the primary. Shocks in the gas might then trigger (1994|) further investigated whether a weak interaction could 



a nuclear starburst as seen in the spectra of Carter et al 



(1988) 



Mergers al so provide th e framework for most model s 
of KDCs (e.g. [Bender 199"c| ). [Pranx fc lUingworth (198^ ) 
noted that a counter-rotating core such as that of IC 1459 
could form in an accretion of a gas rich satellite by gas- 
dynamical evolution and ensuing star formation, or from 
the dynamical disruption of a satellite near the cor e of the 
main gal axy, by simple stellar dynamical processes 
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explain the existence of KDCs, and proposed a mechanism 
whereby a counter-rotating core could be formed by the ret- 
rograde passage of a massive galaxy past a slowly rotating 
elliptical with a pre-existing rapidly rotating central disk. In 
this picture it is the outer regions of the galaxy which are 
anomalous, the direction of rotation is changed or reversed, 
whereas the central disk is relatively unaffected. 

A rathe r different model for the formatio n of shells was 
proposed by Fabian, Nulsen, fc Stewart (198C|) and expanded 



fc Quinn (199[|), testing with JV-body models a scenario first by Bertschinger (1985) and Williams & Christiansen (1986) 



proposed by Kormendy (1984), showed that KDCs includ- 



ing counter-rotating cores could be formed in purely stellar 
dynamical simulations. They studied mergers between el- 
liptical galaxies of unequal mass, and found that the core 
kinematics in the remnant depends mostly upon the orbital 
angular momentum at a late stage of the merger, whereas 
the kinematics of the outer regions are largely the original 
kinematics of the primary. Thus in retrograde encounters 
a counter-rotating core can form. When counter-rotation 



In this picture shells are formed as a result of star formation 
in shocked gas in the interstellar medium of the galaxy, the 
shocks being caused by outflow from the nucleus. With the 
measurement of shell colours, which are not much bluer than 
the main body of the galaxy, the failure to detect either 
ionized or neutral gas associated with the shells except in a 
very few cases, and the discovery of the interleaving of shells 
on opposite sides of the nucleus in man y shell galaxies, this 
picture has falle n into disfavor, though Lowenstein, Fabian, 



forms via this mechanism, both primary and secondary ma- 



terial coun ter-rotate at t he core, and skewed line profiles are 
produced (Balcells 1992). Spiral-Spiral (SS) major mergers 
provide promising mechan isms for the formation of KDCs. 
Hernquist & Barnes (1991) show an SS merger model with 



Nulsen (1987) reconcile this model with the last of these 



observations. 

Motivated by the observational and theoretical argu- 
ments outlined above, we initiated a programme to look for 



KDCs in a sample of shell galaxies drawn from the Malin 



an emb jdded counter-streaming gaseous disk. Although the 



evolution of the gaseous phase is still crude and star forma- 
tion has not been included (hence the predictive power of 
the model is limited), the physics is probably on the right 
track. 



Carter (1983) catalogue, to investigate whether there is a 



common formation mechanism for shells and KDCs. In this 
paper we report the results for NGC 2865. 
NGC 2865 is listed in the 



alogue of shell galaxies. It is classified as E3-I- 



pose a : nechanism for the formation of shells in a weak (hy- 



perbolic) interaction between two massive galaxies. In this 
picture shells are formed in a rotating "thick disk" compo- 
nent of an otherwise non-rotating oblate elliptical, as an one- 
armed ( [ensity wave excited by the passage of another galaxy. 



Malin & Carter (1983) cat- 
RC3 (g^ 



Thjamson fc Wright (1990) and [Thomson (199lD pro- Vaucouleurs et al. 199l|), and E4 in RSA ([Sandage fc Tarn 



mann 1991[ ). Taking mt, = 12.35 from RC3, Ab = 0.27, 
hioo ~ 1 and cz — 2530 kms"'^, its Mb of -20.02 places it on 
the brighter side of the dividing line between low- lum inosity 
and giant ellipticals (Mb^ = —18.99 with hioo = 1; Davies 
et al. 1983[) . CCD photometry reveals surface brightness that 



NGC 2865 3 



is consi stent with a r4 law out to ~ 60 ( JReid, Boisson, fc| Table 1. Instrumental setup for the observations 



Sanson 1994). Its chaotic shells are bluer than the main body 
of the galaxy at ~84-98% confidence level, and cont ribute 
~ 11% - 22% of the total luminosity ( JFort et al. 1986| ). Fort 
et al. noted the colour indices of the shells resemble those of 
an Sb or Sc galaxy. The galaxy appears disturbed in a deep 
R-band image, and in the outer parts the shells are brighter 
t han the main bod y (Fig. hi). The mild "E+A" spectrum 
(Carter et al. 1988) in the nucleus suggests some recent star 
formation. Carter et al. estimated that light from A stars 
contributed 30 ± 10% of the light at the nucleus. NGC 2865 
is detected by IR AS at 100 Jy, with a n estimated dust mass 
of 0.17 X lifMQ ( [Roberts et al. 199l[ ). These properties are 
evidence for a recent merger or accretion that involves a 
gas-rich disk galaxy. 

By performing populatio n synthesis using ga laxies and 



galactic clusters as templates, Bica & AUoin (1987) find that 



a starburst has occurred, and derived an age of 1.2±0.3 Gyr 
for it. 

Schiminovich et al. (1995) observed NGC 2865 with 
the VLA, and found 6 x IQ^K^qqMq of cold gas. Surpris- 
ingly, the gas is found in a broken ring in the outer regions 
of NGC 2865, and is displaced to the outside of the shells 
and the loops. The HI velocity field is smooth, and has the 
same sense and magnitude as the inner stars. Fits of rotation 
curves to the intensity- weighted velocity field yielded a circu- 
lar velocity of 240 ±15 kms~ ^ and an inclination of 65° ± 5° , 
at a position angle of 310° ± 5°, close to the photometric 
major axis. Schiminovich et al. derived a mass-to-light ratio 
M/ Lb of 33 ± 4feioo at large radii, which is about 5 times 
greater than M/ Lb at the core. 

There is a high degree of rotational support for 
NGC 286 5. Stellar kinem atics along 4 axes previously mea- 
sured by Bettoni (1992) shows v/a ~ 1.3, greater than re- 



quired for an isotropic rotator of the same fiattening. Bettoni 
further argued that NGC 2865 must be oblate, because the 
velocity gradients along the 4 position angles are consistent 
with a sinusoidal function of position angle, with the fastest 
rotation along the major axis. Peculiar kinematics were not 

found in the Bett oni data. 

Whilst both Bettoni (1992) and pchiminovich et al 



(1995) point to a merger as responsible for the shell for- 
mation, the physical pictures favored are different. Citing 
the regular stellar velocity field and lack of emission lines, 
Bettoni concludes that his results are consistent with the 



Dupraz & Combes (1986) model in which a small satellite 



is absorbed and disrupted by an almost oblate E3.5 galaxy, 
and that the collisi on has only marginally pert urbed the un- 
derlying potential. Schiminovich et al. (1995), however, fa- 



vor the equal-mass spiral-spiral merger scenario ( Hernquist 
&: Spergel 1992], iHibbard & Mihos 19951) in which the sheUs 



and other fine features are formed out of material raining 
back into the core after the collision, citing as support the 
near spiral Mm/ Lb in the outer regions and the close match- 
ing between the stellar and gas kinematics. They concluded, 
however, that the displacement of the gas outward from the 
stars is not supported by any of the current merger scenar- 
ios. 

In this paper we report on the discovery of a kinemati- 
cally distinct core in the nucleus of NGC 2865. We recover 
the line-of-sight velocity distribution along two axes, and 
find a small, rapidly rotating stellar component with asso- 
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dated ionised gas. We measure the abundance distribution 
of NGC 2865, and find that a young stellar population is 
associated with the KDC. A crude modelling of the dilution 
effects using stellar templates shows that major star forma- 
tion must have stopped at least ~ 0.5 Gyr ago. The KDC 
shows a iron to magnesium index ratio higher than the light 
from the underlying galaxy. By considering the evolution of 
a post-starburst spectrum with time, we find that this may 
be explained as a temperature effect in the atmospheres of 
stars at main-sequence turn-offs between A3 and F4, in- 
stead of a change in abundance ratio. We then constrain the 
star-formation history of the KDC, aided by the spectral 
evolution library of Bruzual & Chariot (1997). We derive 
constraints on physical parameters, and find that the KDC 
was formed either by a metallicity-enhanced starburst, or by 
accretion of a spiral into a giant elliptical with a different 
abundance ratio, during which star formation is truncated. 
The layout of this paper is as follows: in §2 the obser- 
vation and data reduction procedures are described. Stellar 
kinematics derived from the absorption- line profile, and gas 
kinematics derived from the [QUI] emission line are pre- 
sented in §3. In §4 the radial gradients of absorption line 
indices are presented. In §5 & 6 we put initial constraints on 
the age of the starburst by considering the effect of dilution 
by a young stellar population. In §7, the star formation his- 
tory is further constrained using evolutionary stellar popula- 
tion synthesis and considering a "bulge-(-burst" model under 
4 scenarios of star formation. This is followed by a discus- 
sion in §8. In order not to distract from the main scientific 
thrust, the reliability of the methods employed in extract- 
ing the absorption-line profile is investigated by Monte Carlo 
simulations in Appendix H, whilst the dependence on stellar 
templates is investigated in Appendix M. 



2 OBSERVATIONS 

2.1 The data 

Long-slit spectra of NGC 2865 were taken with the RGO 
spectrograph on the Anglo Australian Telescope in April 
1995. The instrumental setup is summarized in Table 111. The 
spectral and spatial resolution are significant irn provements 
over earlier studies of this galaxy (Bettoni 1992). The spec- 
tra span the wavelength range from 4850A to 5610A, cover- 
ing the H/3, Mg & Fe absorption lines as well as the [OIII] 
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emission line at ~ 5007A. With a 2.6 slit, the exposures 
were 3 600 seconds alons the major axis (P. A. 153°) and mi- 



nor axij (P. A. 63°). The seeing was 0.8 FWHM for the 



major axis data and 0.9 for the minor axis data. 

The data were bias subtracted, flat-fielded, cleaned of 
cosmic rays and wavelength calibrated using the standard 
procedures in the AAO FIGARO data reduction package. 
The spectra were rebinned to a logarithmic wavelength scale 
to allow cross- correlation and Gauss-Hermite analysis. 

In § 4.2.2 we find that the central Mg indices are de- 



pressed, but not the Fe. This may arise if there exists a focus- 
ing problem with the spec tralgraph under good seeing con - 
ditions, first mentioned by Davies^ Sadler, fc Peletier fl993) 



and further discussed by Gonzalez (1993 ) and Mehlert et al , 
(1998[K_ This arises as the continuum bandpasses of Mgi and 



Mg2 has a wide wavelength separation, and an out-of- focus 
camera will distribute the fluxes at the end of the central 
spectra away, thus artifically lowering the continuum level 
and in turn the absorption- line index. Such a problem is not 
present in our data. The CCD frames are well and unifor- 
mally focused, with no detectable spatial variation in focus 
when an arc lamp exposure is examined. 

Template spectra of 3 stars HD 80170 (K5III), 
HD 107328 (KOIIIb) and HD 109379 (G5IIb) were taken 
with the same instrumental setup as the galaxy observations. 
After their radial velocities measured by cross-correlation, 
the template spectra were de-redshifted to the reference 
frame of HD 109379, which has a h eliocentric velocity of 
-7.8 kms"^ ([Mayer fc Maurice 1985|). This value is added 



to the measured velocities in subsequent analysis. 



2.2 Gauss-Hermite parameterization of LOSVD 



Following van der Marel fc Franx (1993|), the line-of-sight 



velocity distributions (LOSVDs) L{w) are taken to be Gaus- 
sians with small Hermite deviations: 



L{w) — J- 



where 



q(xu) 



t(ui) 



1+^/ln Hni^ 



V 



1 _™i 

e 2 



2-K 



(1) 

(2) 
(3) 



and 7, V, a are the line-strength, mean velocity, velocity 
dispersion respectively. H„ is a Hermite polynomial of or- 
der n and with a corresponding amplitude hn- Expressed 
explicitly, H3 and H4 are: 






3w) 



(4) 
(5) 



Thus the LOSVD is parameterized by the variables (7, V, 
a, /13, ..., /ijv), and its shape is a perfect Gaussian when 
the hn are zero. The Gauss-Hermite moments h^ and /14 
are meaningful parameters as they are related to the skew- 
ness and the kurtosis of the distribution respectively, and 
are useful indicators of unusual core kinematics. A large /13 
with opposite sign to the mean velocity is typical of rapidly- 



rotating systems, whilst /14 is relate d to the velocity disper- 
sion anisotropy in the galaxy core (van der Marel fc Franx 



199J). Although the higher moments are more difficult to 
interpret, they are still useful for identifying structures in 
the residuals of the model fits, and furthermore a large sig- 
nal in them indicates the breakdown of the near-Gaussian 
approximation, therefore parameters up to he are measured 
and plotted. 

In the rest of this section, the method of extracting the 
LOSVD is described. Readers not interested in this should 
go straight to § Bjl. To recover the LOSVD, we use the pro- 
gram kinematics, written and kindly made available to us 
by H.- W. Rix. This program impleme nts the parameteriza- 
tion of vanderMarel&FVanxfl993 ) into t he pixel space 
fitting routin es of Rix fc White (1992 ) ; see also yan der Marel 
et al. (1994). First a stellar template is convolved with an 
initial guess of L{w) to give a model galaxy spectrum. The 
parameters (7, V, a) are recovered when the noise weighted 
X^ of the model and the galaxy spectrum are minimized. 
Each hn is estimated by keeping all lower moments fixed in 
the x^ fitting, and this process repeated until n = 6. r] The 
errorbar of each parameter is taken as the 63% confidence 
region, assuming all other parameters are fixed. 

The wavelength range 4910-5503 A is well known to 
give template mismatch problems. To minimize this, an op- 
timal stellar template is employed, taken as the linear combi- 
nation of the 3 templates that minimizes the x^ of a Gaus- 
sian profile, after first estimating (7, V, a). The optimal 
template, derived separately for each point along the slit, is 
used to recover the LOSVD, modelled up to he. The galaxy 
nucleus is located by a Gaussian fit to the light profile of the 
central few pixels. Due to the presence of a foreground star, 
no results are presented beyond —6" of the nucleus along 
the minor axis. 

Experimentation shows that summation to continuum 
S/N > 40 per pixel gives satisfactory results. The sensi- 
tivity of the measured parameters to noise and template 
mismatch are investigated in Appendices H fc H. The con- 
clusions from Monte Carlo simulations in Appendix H are 
the following: Hermite moments up to fee can be recovered 
to within 0.01-0.02 up to 3 from the nucleus. A small but 
significant displacement in a parameter is introduced if the 
Hermite moment two orders up is large. The parameter least 
sensitive to noise is v, whilst a can be affected by S/N and 
by power in /14, with a signal of 0.05 in /14 causing a shift of 
~ -)-4kms~'^. There are no problems associated with imper- 
fect sky subtraction, even when 12 spectra are combined to 
achieve the desired S/N. Therefore along the major axis the 
LOSVD is modelled up to he for \r\ < 3 and up to hs for 
|r| < 6 , and to a smaller extent along the minor axis due 
to the lower S/N. 

The conclusions from Appendix H are the following: 
Template mismatch cannot explain the radial profile of the 
parameters (apart from possibly h^). The maximum offsets 
are of order 6kms-^ in v fc cr, and 0.026, 0.026, 0.020 fc 



* In this paper we give {v, 
L(w) according to equation 



) as those defining the best fitting 
, which are slightly different to (V , 



a) defined as moments of V and V ^ of a LOSVD, as pointed out 
by van der Marel fc Franx (1993|), who also give expressions of 



the formulae for conversion. 
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0.018 in /13, hi, h^ & he respectively. The observed zero- 
point offsets in ^3, /15 & /is may be attributed to template 
mismatch, but the positive offset in ht is marginally signif- 
icant even in the extreme case. Thus we think the positive 
hi reported below has a physical origin. 



2.2.1 Regions excluded from direct fitting 

Emission lines can potentially introduce systematics in the 
absorption line profile analysis, and further affects the mea- 
surement of absorption-line indices if they fell within one 
of the bandpasses. Hidden [NI] emission at 5199A can arti- 



ficially lower the Mg h equivalent width by 0.4-2 A (Goud- 
frooij fc Emsellem 1996|), but given that no significant resid- 



uals can be seen at 5199A from the X^ fits, it is probably neg- 
ligible. From the x^ fits we find some extremely weak resid- 
uals within the continuum bandpasses of Mg b and Fe527o 
coincide with the location expected of some [Fe] or [CrI] 
emission lines. Whilst these residuals may be caused by tem- 
plate mismatch, nevertheless the regions 4951. 0-4955. 9A, 
5141.27-5158.lA & 5519.6-5533.8A (in the galaxy's rest 
frame) are excluded from the line-profile analysis. The H/3 
absorption line, [OIII] (5002.7-5014.lA) and [NI] (5196.4- 
5201.5 A) emission lines are also excluded, as well as oc- 
casional regions affected by cosmic rays. Experimentation 
shows that excluding a few pixels from the x^ fits does not 
affect the recovery of the kinematic parameters because in- 
formation is derived from the entire spectrum. 



3 CORE KINEMATICS OF NGC 2865 

3.1 Kinematics of stars derived from absorption 
line profiles 



geometry along the 2 axes. For the major axis data each pixel 
samples 0.7 (0.09/i^qo kpc) along the major axis and 2.6 
(0.33/i^gQ kpc) along the minor axis, and vice versa for the 
minor axis data. If the kinematic subcomponent is a rapidly 
rotating disk aligned along the major axis and its axis in- 
clined at an angle i with respect to the line-of-sight, then 
the major axis pixels sample mainly the coherent rotation 
of the stars, with the dispersion mainly from the a^cos i com- 
ponent. In comparison, the minor axis pixels are 3.7 times 
wider along the disk and 0.27 times narrower perpendicular 
to it, hence a major contribution to the nuclear a is from 
the velocity gradient over 0.33/i]~qq kpc along the disk, and 
the contribution from a^cosi is smaller. Thus if the disk is 
cold and rotating rapidly it is possible to measure a different 
nuclear o along the major and minor axes. 

Along both axes, the line profile is significantly 
Lorentzian with hi up to 0.05. There is a small offset of 
~ 0.02 in /15 and he along both axes, albeit these 2 param- 
eters are more sensitive to noise and template mismatch. 
Appendix f\ shows that a significant signal in ha of ~ 0.05 
may artificially enhance /14 by ~ 0.02, whereas an /14 of 0.05 
may artificially shift a upwards by about 4kms~^. There- 
fore if the observed /14 and /ig have a physical origin, the 
true hi and a may be lower than measured by about 0.01 
and 4kms~'^ respectively. 

The central positive /14 is one of the outstanding prop- 
erties of the nuclear kinematics of NGC 2865. Our error 
estimates (Appendix A & B) indicate that the signature is 
real, unless there is an extreme template mismatch, unlikely 
given the good match in metallicity between the templates 
and the galaxy. Positive /14 is normally attributed to pre- 
dominance of radial orbits. However, radial orbits yield an 
outwardly decreasing o profile for a constant M / L, which is 
not seen at the nuclear region. We find a similarl y positive 
hi signature in ano1: her shell lenticular NGC 474 ( Hau, Bal 



The re sults of the absorption line profile analvsis are pre- 



cells, & Carter 1996). Su ch positive hj values are rare in th e 



sented in Fig. |. Along the major-axis the rotation curve nuclei of ellipticals (e.g. [Bender, Saglia, fc Gerhard 1994D 

is symmetric, and is near solid-body for jrj Si 3 , reaching 
~ ±63kms-i 

inward from > 

-1 



at ~ 12 . The velocity dispersion increases 
' 186kms~^ at \r\ « 12 to a maximum of 
200kms~^ at |r| ~ 6 . The h^, is consistent with ~ 0.02 
beyond ~ 4 . Within 3 of the nucleus, the rotation is en- 
hanced relative to the general trend, with associated a that 
is depressed by ~ 23kms~^, and point-symmetric /13 with 
magnitude up to 0.05. The latter indicates the line profiles 
are skewed in the opposite sense to the mean velocity. These 
results suggest that NGC 2865 harbours a cold, rotating nu- 
clear stellar component that is co-rotating with the under- 
lying elliptical, giving rise to the steeper velocity gradient 
and point-symmetric /13 in the inner 3 . This component 
has high surface brightness and the light contributed by it 
lowers the central velocity dispersion. Along the minor-axis 
the rotation has an upper limit of 15kms~^, and there is 
no significant signal in h^ except a zero-point of -1-0.01. The 
lack of rotation alon g the minor ax is is consistent with the 
kinematic profiles of Bettoni (1992| ), and suggests that the 
intrinsic shape of NGC 2865 is close to oblate. 

The minor axis a within 3 of the centre are ~ 9 kms~^ 
higher than those in the same region along the major axis. 
The region of lowered a is slightly more extended and shal- 
lower, and the depression of ~ 18kms~^ is smaller than the 
major axis. This difference could be due to the different pixel 



and could contain information on the type of merger which 
produced the shells. 

The central velocity dispersion is about 200 kms^^ with- 
out the KDC. Taking the rotation ve locity of ~ 24 0kms~^ 
at the outermost points of Fig. lb in (Bettoni 1992), we ar- 
rive at a v/a ~ 1.2, slightly smaller than the value of 1.3 
estimated by Bettoni 1992. This difference can be attributed 
to the central dispersion value of 180 kms~^ adopted by Bet- 
toni, who did not detect the KDC. We mea sure a mean e l- 
lipticity of ~ 0.26 for r > 30" . From Fig. 1 of[Binney (1978|), 



we find that the v/a of NGC 2865 is twice that expected for 
an oblate isotropic spheroid of the same flattening. 

3.2 Kinematics of ionized gas 

In this section we present the line-of-sight velocity distri- 
bution of the ionized gas inferred by the profile of the weak 
[OIII] emission line, obtained directly from the residual map 
of the Gauss-Hermite analysis. We find this method very 
successful and more objective than fitting the baseline by 
hand. The results for the 2 axes are plotted in Fig. 0, with 
the stellar rotation curve from the previous section over- 
plotted. Because the [OIII] emission is weak, we feel that 
Fig. p^ allows only a qualitative description. At the nucleus, 
light from the underlying galaxy contributes up to ~ 4700 
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Figure 2 

S/Af > 40 



The V, (7, /i3, /i4, /15 & hg along the major and minor axes of NGC 2865. The results with solid dots are from spectra with 
per pixel, whilst the point plotted with an open symbol has S/N ^ 30. The Icr errorbars are also plotted. 
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,s kinematics from the [OIII] omission line at 5007A along major and minor axes of NGC 2865. The stellar velocity 
are over-plotted in white. No data is plotted for r < —6 on the minor axis due to light contamination by a foreground 



counts per pixel, therefore the S/N per pixel is about 3, 
explaining the noisier appearance near the nucleus. 

Along the major axis, the gas co-rotates with the stars. 
At the location where the kinematically distinct compo- 
nent is observed in the stellar component, the [OIII] line 
appears to be split into a parallelogram shape with a peak- 
to-peak velocity of ~ 260kms~^. Along the minor axis there 
is no observable rotation. Within 1 of the nucleus, the 
[OIII] line appears to have 2 peaks — at ~ 2570 kms~^ and 
~ 2780 kms~^, but we note that the data along the minor 
axis have a lower S/N. As a rough estimate, assumi ng th e 
gas is in circular orbits which are inclined at 45° (§ 4.2.1), 
and taking a projected rotational velocity of 130 kms^^ 
along the major axis, the total mass within 1.9 (0.24 /i^gg 
kpc) of the nucleus is estimated to be 1.9 x 10^h~^MQ. 

Despite the low S/N, the major and minor axes [OIII] 
brightness are well fit by power-law profiles, and not an ex- 
ponential profile. The power-law indices for the major and 
minor axes profiles are 0.651 and 0.891 respectively, indicat- 
ing the light profile is steeper along the minor axis, and that 
the flattening along the major axis increases with distance 
from the nucleus. 

The results from this section show that the stars and gas 
rotate with approximately the same magnitude and direc- 
tion about the minor axis. The flattening and lack of rotation 



along the minor axis suggests that the gas is a co-rotating 
disk asiJoclaLed with the nuclear kiuemalic sub-couipoueuL. 



4 ABSORPTION LINE INDICES 

4.1 Procedures 

In this section the radial profiles of the absorption line in- 
dices are presented. Line indices are extracted from individ- 



ual, de-redshifted spectra following the recipe of Faberetal 



(1985), acc ording to the bandpass definitions of Worthey 
et al. (1994| ) adopted in the Bruzual & Chariot models (§0) 



whilst the molecular bands Mgi & Mg2 are measured in 
magnitudes. The indices are corrected to a zero dispersion 
system by applying correction factors, estimated by broad- 
ening the 3 stellar templates to dispersions ranging from 
160kms~^ to 210kms~^. The correction factor for each a 
is taken as the mean of the correction estimated for the 3 
stars. The typical correction factor at ct = 200 kms^^ for 
Mgi, Mg2, Mg b, H/3, Fe527o and Fesass, are 0.976. 0.994, 
0.873, 0.968, 0.886 & 0.865 respectively (where the measured 
line indices should be divided by these values). The spread 
in these values ac are 0.006, 0.005, 0.009, 0.028, 0.012 and 
0.009 respectively. P] The quoted errors aj in the measured 
line indices take into account both the error due to photon 
shot-noise and ac- The indices are not degraded to the Lick 
spectral resolution of ~ 8A FWHM, and because our data 
are not flux calibrated, and no stars with Lick indices are 
observed, there is some uncertainty in the zero points. This 
does not affect the interpretation of our results; firstly be- 
cause the Mgi and Mg2 indices (which show radial variation) 
are extremely insensitive to line-broadening, and secondly 
because the difference in each index (e.g. AH/3, AFessss & 
A<Mg> = (Mgi + Mg2)/2) is the free parameter in the 
models. Furthermore, within 8 of the nucleus along the ma- 
jor axis, the average Mg2 of 0.192 ± 0.014 mag agre es with 
that measured by the Lick group (0.208 ± 0.02 mag; Trager 

(1997| )). I The average Fe5270 and Fessss of 2.75 ± 0.24 A 



Passbands for the various indices and continua are given 
in table H. The atomic absorption lines II/3, Mg b, Fe527o 
and Fe5335 are measured as equivalent widths in Angstroms, 



and 2.77±0.17 A also agree within error with the Lick group 
values (3.12±0.4 A and 2.82±0.4 A). Thus we are confident 
that our measurements are close to the Lick scale, and the 



T The dispersion correction factors for H/3 arc estimated by re- 
defining the central bandpass to 4861-4877A and using only the 
red continuum level, because at zero redshift the blue continuum 
bandpass of the template stars is shifted out of the observing win- 
dow. This method accounts for a larger ac compared to other line 
indices, but is not a problem as the error in H/3 is dominated by 
photon shot noise. 

+ A proper comparison would require convolving the window 
function adopted by the Lick group with our data. 
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Figure 4. The radial profiles of the absorption line indiees Fe527o, Fe5335, H/3, Mgi, Mg2 and Mg h along the major and minor axis of 
NGC 2865. 
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Table 2. The bandpass definitions of measured absorption line indices, according to Worthey et al. (1994). 



Feature 


Type 


H/3 


Atomic absorption line 


Mgi 


Molecular band 


Mg2 


Molecular band 


Mgb 


Atomic absorption line 


Fe5270 


Atomic absorption line 


Fe5335 


Atomic absorption line 



Central Bandpass (A) Continuum Bandpasses (A) 



4847.875-4876.625 
5069.125-5134.125 
5154.125-5196.625 
5160.125-5192.625 
5245.650-5285.650 
5312.125-5352.125 



4827.875-4847.875, 
4895.125-4957.625, 
4895.125-4957.625, 
5142.625-5161.375, 
5233.150-5248.150, 
5304.625-5315.875, 



4876.625-4891.625 
5301.125-5366.125 
5301.125-5366.125 
5191.375-5206.375 
5285.650-5318.150 
5353.375-5363.375 



zero-point offsets are the same order of magnitude as the 
smaUest errorbars. 



4.2 Radial profiles of absorption-line strengths 

4.2.1 HP 

The radial profiles for the Fe527o, Fessas, H/3 Mgi, Mg2 & 
Mg b indices are presented in Fig. y. The H/3 is ~ 2 A at large 
radii, and increases gradually towards the nucleus, where 
it reaches ~ 3.5 A. H/3 absorption is strongest in A- type 
stars and a large equivalent width indicates the presence of 
a young stellar population. The major axis H/3 profile in- 
dicates that the fraction of the young population along the 
line-of-sight gradually increases towards the nucleus, with a 
possible further enhancement where o is depressed, suggest- 
ing that the young population is associated with the kine- 
matic subcomponent. Along the minor axis the H/3 profile is 
steeper and occupies a smaller region, indicating the extent 
of the young population is a factor ~ 0.7 smaller along the 
minor axis. If the young population is confined to a disk, 
its angle of inclination would be about 45° with respect to 
the line-of-sight, and curiously would have approximately 
the same orientation as the HI ring in the outer regions 



( Schiminovich et al. 1995 ). The reason why the possible dis- 
continuity of the H/3 index observed along the major axis is 
not apparent along the minor axis may be the caused by a 
smaller extent of the young population along this direction, 
and/or a larger seeing. 

The presence of weak [OIII] emission in our spectra sug- 
gests there may be some hidden H/3 emission too, which 
could lower the measured equivalent width, therefore the 
measured H/3 at the nucleus is probably a lower limit. 



Similar profiles are observed along the minor axis, but with 
a noisier appearance due to lower S/N . 

The absorption line indices of the bulge could be low- 
ered if the continuum level is raised by light of very young 
stars. The enhanced H/3 and lowered Mg indices are, in the 
first instance, consistent with the picture in which the nu- 
clear light originates from a young stellar population associ- 
ated with the kinematically cold subcomponent, and from an 
old population associated with the kinematically hot bulge. 

The observed Mg2 indices are more than 0.08 mag lower 
than that expected from the Mg-a relationsh ip, which pre - 
dicts Mg2 ~ 0.294 mag for a = 200 kms"^ ( |Davies 1996| ). 
This difi'erence is too large to be accountable by the pos- 
sible difi'erence in zero-points. Together with the high H/3 
(~ 2 ± O.4A) of the "bulge", this may suggest that some 
young stars are not merely confined to the kinematic sub- 
component, but is also mixed into the dynamically hot bulge 
population. 



4.2.3 Fe 

The Fe527o and Fesass profiles are quite different to those 
of Mg, with no observable gradient in either index. This 
result is rather surprising given that there are pronounced 
depressions in the Mg indices. If the depressed Mg indices 
at the nucleus are due to dilution from a young population 
with featureless, blue spectra, then the Fe is expected to 
be depressed too. Although Fe527o could be enhanced by 
hidden emission in the continuum bandpass, no significant 
depression is seen in Fessss , therefore the flat profiles cannot 
be attributed to hidden emission. In the next section we shall 
quantify the amount of depression expected in the Fe indices 
due to light from a young stellar population. 



4.2.2 Mg 

The Mg indices show a gradual increase towards the nucleus, 
but they are lowered sharply at the location of the kinemati- 
cally distinct core. Along the major axis Mg2 increases from 
~ 0.18 mag at r « 16 to ~ 0.205 mag at r ~ 4 , but then 
lowers to ~ 0.153 at the nucleus. The Mgi rises from ~ 0.06 
mag at r « 16 to ~ 0.09 mag at r ~ 2 , but then low- 
ers to ~ 0.037 at the nucleus. The same trend can be seen 
in the Mg 6 equivalent width, albeit with larger errorbars. 
The Mg b also has a different shape to the Mgi and Mg2 at 
the nucleus, and the depression is shallower and less sharp. 
As discussed earlier in 



2.2.1 



Mg b is more prone to hidden 
emission at the nucleus and is probably an upper-limit there. 



5 CAN THE NUCLEAR LINE INDEX 

GRADIENTS BE EXPLAINED BY LIGHT 
DILUTION BY 05-A3 STARS? 

We want to investigate whether light dilution by young O, 
B or A type stars alone are able explain the observed line- 
strength gradients at the nucleus of NGC 2865. For this pur- 
pose, the spectrum at —3.37 on the major axis of NGC 2865 
is extracted, and varying proportions of spectra between 
type 05 and A3 of class I, HI and V taken from a library 
of fiux-calibrated, solar-metallicity spectra are added to it 
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0,15 0.14 0.13 0.12 0.11 0.1 0.09 

<Mg> (Mag) 

Figure 5. Effect of light dilution on the Fesaas line by adding 
spectra of O, B & A type stars selected from the spectral library 
of Jacoby, Hunter & Christian (1984). The data points along the 
major axis of NGC 2865, with |r| < 4 , are plotted. The spectrum 
at r = —3.37 is extracted and various amount of stellar spectra 
between 05 & A3 are added to it. The arrows denote the trend 
in which the <Mg>=(Mgi+Mg2)/2 and Fessas follows as the 
galaxy spectrum is diluted more and more by light from the early- 
type stars. 



( Jacob^fT^imil^i', &i Clii'liJLlaii_1984 ). p] Tliti lltie-lMl(;eiJ 



^ , ^ ^ ni'y 

then measured, and in Fig. |5| the trend which Fessas and 
<Mg>=(Mgi+Mg2)/2 follow as a result of light dilution 
are plotted over the actual data points at \r\ < 4 along 

the major axis. P] The equivalent width of the Fesaas line is 
expected to be lowered by at least 0.55 A for a change of 
~ —0.052 mag in <Mg>. This depression is not observed, 
and it is clear that the data points do not follow the trend 
expected from dilution. Therefore the primary component 
of the nuclear light in the core of NGC 2865 does not orig- 
inate from 05-A3 stars, and the gradient of the spectral 
continuum cannot account for the different behavior of the 
Mg and Fe. In order to explain the observed line-strength 
gradients, one must appeal to a different evolution of in- 
dividual metallic lines as the starburst ages — this shall be 
investigated in the next section. Results in this section also 
demonstrate that major star formation must have ceased for 
at least o(5 x 10*) yr, the lifetime of a Main Sequence A3 
star. 



S Although our spectra are not flux calibrated, this will not 
change the conclusions, firstly because the gradients of the lines in 
Fig. Hare independent of the zero point, and secondly because the 
separation between Fe and Mg absorption lines is small compared 
with the wavelength range covered. 

'I We chose to model <Mg> to boost signal and to reduce the 
number of parameters to model. Fesaas is chosen because Fe527o 
could be affected by hidden emission. As it turns out, similar 
results are obtained when <Fe>=(Fe527o+Fe5aa5)/2 is modelled 
instead. 



6 HOW DO THE MG AND FE ABSORPTION 
LINES VARY WITH THE MAIN SEQUENCE 
TURNOFF? 

During the first 8 Gyr after a starburst, the Main Sequence 
(MS) light contribution to the V-band luminosity oscillates 
with t ime between 30 and 75% (assuming Salpeter IMF; 
Fig. 8 [Chariot fc Bruzual 1991 ). Furthermore, during this 
time the spectral energy distribution of the MS will approx- 
imate that of the stars at the turn-off, because their bright- 
ness outweighs the number of stars at later types. How do 
the Mg and Fe indices behave as a function of spectral type? 
Line indices are measured for a sampl e of standard solar - 
metallicity stars with spectral class V ( Jacoby_et al. 1984 ) 
and plotted against the spectral type in Fig. H. Spectra of 
MS stars earlier than A3 do not have a significant signal in 
<Mg> and <Fe>. Between approximately A3 and F4, the 
<Fe> is significant, whilst <Mg> is still negligible. Only 
after about F4 are both <Mg> and <Fe> significant. 

The ages of an instantaneous burst with MS turn-off 
at A3 and F4 are estimated in the following way: Taking 
(B-V) of solar metalli city A3V and F4V star s to be 0.09. 
and 0.39 respectively (Mihalas & Binney 1981), their corre- 
sponding e ffective temperatures are 8760 K and 6700 K re- 
spectively ([Mihalas fc Binney 1981], [Alonso et al. 199e|). By 



linearly inte rpolating the theoretical isochrones of 



BerteUi 



et al. (1994), we arrive at age estimates of 0.5 Gyr and 1.7 



Gyr for MS turn-offs at A3 and F4 respectively. This result 
is consistent with Bica & AUoin's (1987) age determination 
using star cluster templates. 

In summary, about 0.5-1.7 Gyr after an instantaneous 
burst, the Fe absorption lines in the burst population may 
be stronger than that of the Mg. This is due to a temper- 
ature effect in the atmospheres of stars at the MS turnoff. 
When the spectrum of this population is added to that of 
the old population, the total <Mg> will be lowered more 
than the <Fe>, and the inferred [Mg/Fe] abundance ratio 
appears to be modified. This difference may be enhanced by 
increasing the light contribution from the burst population 
(i.e. by increasing its strength), by increasing the overall 
metallicity of the burst, or by increasing its Fe relative to 
Mg. Indeed, population synthesis in the next section shows 
that an enhancement in metallicity or Fe abundance is re- 
quired to explain the nuclear line indices of NGC 2865. Note 
that this effect only applies if the age spread of the newly 
formed stars is smaller than the MS age difference between 
A3 & F4 (i.e. << 1.2 Gyr). We also note that although 
the difference in the Mg and Fe indices can be attributed 
to the temperature effects of the stars at MS turnoff, if the 
metallicity is high enough, after about 1 Gyr the Mg and 
Fe indices in the burst population are comparable to those 
of the old population, thus kinematics measurements can 
be performed (§ 3.1). In the next section, the observed in- 



dices are compared with predictions by evolutionary stellar 
population synthesis. 



7 CONSTRAINTS ON THE RECENT STAR 
FORMATION HISTORY BY STELLAR 
POPULATION SYNTHESIS 

We shall constrain the recent star formation history of 
NGC 2865 using population synthesis, by considering a 
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Table 3. A summary of the parameters adopted in the star formation models. The first 6 columns are, respectively: model name, the 
start and stop time of the secondary episode of star formation in Gyr, the metallicity of the bulge and the nuclear subcomponent, and the 
[Fe] overabundance (in dex) of the nuclear component relative to solar. The remaining 2 columns summarize the model parameters which 
satisfy the observational constraints. They are the time elapsed since cessation of the secondary episode of star formation At = inom —tstop, 
and the "burst fraction" /;, along the line-of sight at the nucleus. 



Model 


tstart 


istop 


^bulgei'Z'Q) 


Znuci'Z'Q) 


A [Fe] (dex) 


At (Gyr) 


fb 


SBl 


7 


7.1 


1 


1 





_ 


_ 


SB2 


18 


18.1 


0.4 


0.4 





- 


- 


SB3 


7 


7.1 


1 


2.5 





0.7 
0.8 
0.9 
1.0 


0.33 
0.50 

0.67-0.91 
0.91 


SB4 


18 


18.1 


0.4 


1 





1.0 
1.2 
1.3 
1.4 
1.5 
1.6 
1.7 


0.20 
0.33 
0.33 

0.50 
0.50 

0.67-0.91 
0.91 


SB5 


18 


18.1 


0.4 


2.5 





0.4 
0.5 
0.6 
0.7 
0.8 
0.9 


0.13 

0.20-0.33 

0.20-0.91 

0.33-0.91 

0.50-0.91 

0.91 


SB6 


7 


7.1 


1 


1 


0.08 


- 


- 


SB7 


7 


7.2 


1 


2.5 





0.6 
0.7 
0.8 
0.9 
1.0 


0.20 

0.33-0.50 

0.50-0.91 

0.67-0.91 

0.91 


SB8 


7 


7.4 


1 


2.5 





0.5 
0.6 
0.7 
0.8 
0.9 


0.20 

0.33 

0.50-0.67 

0.67-0.91 

0.91 


TSFl 


1 


7 


1 


1 





- 


- 


TSF2 


1 


7 


1 


1 


0.04 


- 


- 


TSF3 


1 


7 


1 


1 


0.08 


0.1 
0.2 
0.3 
0.4 


0.50 

0.67-0.80 
0.80-0.99 
0.99 


CSFl 


1 


- 


1 


1 





- 


- 


CSF2 


1 


- 


0.4 


0.4 





- 


- 


RCSFl 


7 


- 


1 


1 





- 


- 


RCSF2 


7 


- 


1 


2.5 





- 


- 



"bulge+burst" model in which nuclear light originates from 
an old bulge population plus a younger population associ- 
ated with the kinematically distinct component .[U Further- 
more, the star formation histories are simplified by consid- 
ering 4 basic shapes, illustrated in Fig. Q. 

Due to degeneracy, the details of star formation cannot 
be distinguished after a few Gyr. So for the purpose of this 
work, the bulge stars are formed between & 1 Gyr after 
the Big Bang with a constant star formation rate SFR. Ini- 



tially, the bulge H/9 is strong due to the dominance of young, 
bright early-type stars, whilst the metallic line-strengths are 
weak. As the bulge ages, its H/3 falls whilst its metal line- 
indices rise and approach some asymptotic values similar to 
those observed at present. The introduction of a secondary 
episode of star formation (the "burst") immediately lowers 
the metallic line indices of the combined nuclear light and 
raises its H/3, due the presence of early-type stars. Whilst 
star formation is ongoing, H/3 remains high rj, whilst the 
metallic line indices remain lowered. After star formation 
stops, the line-indices will again approach some asymptotic 



II For convenience, throughout this work we shall refer to the 
separate, younger population as the "burst", even though it may 
have prolonged star formation. 



** Ignoring the presence of emission lines when the star forma- 
tion rate is high 
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Figure 6. The <Mg> = (Mgi +Mg2)/2 and <Fe> = (Fe5270 + 
Fe5335)/2 line-strengths as a function of spectral type for Main 
Sequence stars selected from the library of Jacoby et al. (1984). 
Where more than one star of the same spectral type is available, 
an average and its 1 a errorbar is given. 



values. The H/3 of the combined nuclear light approaches 
that of the bulge from above. The asymptotic values of the 
metalhc line indices depend on the metallicity Znuc of the 
new stars. In the absence of self-enrichment, they would be 
higher than the bulge values if Znuc is higher than the bulge 
metallicity Zi,uige, and vice versa. 

The absolute values of our line-indices cannot be mod- 
elled due to uncertainty in their zero-points. Therefore we 
model the line-index differences between the combined nu- 
clear light (i.e. bulge-|-burst) and that of the bulge. Fur- 
thermore we concentrate on modelling AH/3, A<Mg> & 
AFe5335 , because Mg b an d Fe5 27o could be affected by weak 
emission, as discussed in 



2.2.1 



We adopt A<Mg> between 
the nucleus and the bulge to be A<Mg> = -0.053 ± 0.012 
mag, and AFe5335 = 0.00 ± 0.35A. Because of possible hid- 
den H/3 emission, indicated by weak [OHI] emission, AH/3 is 
taken to be AH/3 > 1.4A. Star-formation rates which satisfy 
these observational constraints are sought. 

Star formation rates (SFR) of the 4 basic shapes il- 
lustrated in Figra are investigated, with the adopted pa- 
rameters summarized in Table H. In models SB1-SB9 we 
investigate whether the nuclear subcomponent is formed in 
a recent starburst which r nay be triggered by, for e xample, 
a spiral-elliptical merger ( Hernquist fc Weil 1992 ). These 



models can also apply to spiral-spiral mergers as long as 
the stars are well segregated into two distinct populations, 
that of the older and dynamically hot "bulge" consisting of 



Short Burst 



tstop 



time 



SFR 



SFR 




time 



A "■ Start 



Continuous star formation 



-H* time 



SFR' 



Recent continuous star formation 



time 




Figure 7. The 4 different star formation scenarios considered. 
The abscissa is the age of the galaxy and the ordinate is the star 
formation rate. The star formation rate of the galaxy bulge is 
shaded in grey and that of the nuclear component is shaded in 
black. 



stars formed before the collision, and the young compact 
nuclear "burst" population formed from gas during the col- 
lision. We note in this case the bulge stars on average would 
appear to be younger than the Universe. Next, in models 
TSF1-TSF3, we investigate whether the nuclear stars were 
accreted from a star-forming galaxy in an event which also 
truncated the star formation in this population; the feasibil- 
ity of this model is discussed in § pi Next, in model CSFl, we 
investigate whether the nuclear subcomponent is primordial 
and formed by constant low-level star formation over a Hub- 
ble time up to the present. This model provides a simple ap- 
proximation to small embedded disks in a "disky" elliptical, 
proposed to be a continu ation of the morphologic al sequence 
of spirals into ellipticals (Scorza & Bender 1990). Finally in 
models RCSFl & RCSF2, we investigate whether the nu- 



clear component was formed recently by continuous and on- 
going infall of cold gas that could be returned tidal material 
after a spiral-spiral merger (Hibbard & Mihos 1995). The 
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range of model parameters which satisfy the observational 
constraints are summarized in the last 2 columns of Table pi 
The Galaxy Isoch rone Synthesis Spectral Ev olution Li- 



brary (GISSEL96) of Bruzual fc Chariot (1997) is used as 



a diagnostic tool. It follows the evolution of a stellar popu- 
lation along the evolutionary tracks of the Padova School, 
up to an age of 20 Gyr, with non-solar metallicities mod- 
elled by Kurucz's (1995 version) model atmospheres. The 
library we use is restricted to metallicities of 0.4, 1 & 2.5 Z0 
([Fe/H]=-0.330, 0.093 & 0.560 respectively) and a Scalo 
IMF. Although in reality the IMF depends on the details 
of star-formation, there is still little information in this sub- 
ject area, therefore we use a Scalo IMF and discuss the likely 
effect of changing it later in § ra. 

To match the bulge Mg and Fe line-indices to the ob- 
served levels, a bulge age of about 8 Gyr is needed for 
ZbuigB = 1 Zq, and about 19 Gyr for Zbuige = O.4Z0. Be- 
cause of uncertainties in the age of the Universe, bulges with 
metallicity of l.OZ© & O.4Z0 are modelled. 



7.1 Recent, short burst of star formation 
(SB1-SB8) 

Could the nuclear subcomponent be formed in a recent star- 
burst? We form the nuclear subcomponent at tatart in a star- 
burst which lasted 0.1 Gyr. The metallicity of the burst pop- 
ulation can be the same or higher than that of the bulge. A 
subcomponent with Znuc < Zbuige is not modelled because 
large A<Mg> & AFessas are achievable with low SFR and 
gives small AH/3. The starburst strength is parameterized 
by burst fraction /;,, defined as the ratio of burst mass to 
total mass of the system along the line-of-sight. 



7.1.1 Starbursts with the same metallicity as the bulge 
population (SBl, SB2) 

In Fig. 1^ the line-indices for model SBl where Znuc = 
Zbuige = 1 Zq are presented. Shaded are also the observa- 
tional constraints. The observed AH/3 & A<Mg> are reach- 
able within 1.0 Gyr after the starburst, but the predicted 
|AFe5335| ^ O.5A is ruled out by observations at 92% con- 
fidence level. We find model SB2, where the metallicity is 
0.4 Zq, is also ruled out by observations, and the results are 
not plotted. 



7.1.2 Starbursts with metallicity higher than the bulge 
population (SB3-SB5) 

The above models are repeated with Znuc > Zbuige. The 
line-indices are presented in Fig. 0, where the top row of 
diagrams corresponds to Zbuige = 1 Z0 & Znuc = 2.5 Z©, 
the middle row corresponds to Zbuige ~ 0.4 Z0 & Znuc = 
1 Zq, and the bottom row corresponds to Zbuige ~ 0.4 Zq 
& Znuc = 2.5 Zq. The best-fit parameters are summarised 
in Table |[ 

Since Znuc > Zbuige, the metal line-indices after the 
starburst rise above the asymptotic values of the bulge. The 
Fe5335 rises faster than <Mg>, and crosses over the bulge 
value earlier, at 0.9-2.5 Gyr after the starburst, confirming 
our findings in § |[ If Zbuige = IZq and Znuc = 2.5 Zq, 
solutions exist some 0.7-1.0 Gyr after the starburst, for 



0.33 < fb< 0.91. If At = tnou, - tstop = 0.7 Gyr, a fb = 0.33 
is required. If At = 0.9 Gyr, 0.67 < fb < 0.91 is required. If 
At = 1.0 Gyr, a high /;, — 0.91 is required. 

If Zbuige = 0.4 Zq, SolutloUS exist for both Zburst = 

1 Zq & 2.5 Zq, with the latter satisfying the constraints 
better. If Zburst = 1 Zq, the observational constraints can 
be satisfied if a burst with 0.2 < /& < 0.91 stopped between 
1.0 & 1.7 Gyr ago (see Table H). If Zburst is raised to 2.5 Zq, 
a more recent burst 0.4-0.9 Gyr ago with 0.13 < /[, < 0.91 
is needed. If At = 0.4 Gyr, /,, = 0.13 is required. If At = 0.6 
Gyr, 0.20 < fb < 0.91 is required, whilst if At = 0.8 Gyr, 
0.50 < fb < 0.91 is required. Finally, if At = 0.9 Gyr, a 
rather high fb of 0.91 is required; interestingly at this time 
the nuclear Fe5335 is ~ O.lA higher than that of the bulge. 
In summary, solutions exist if the nuclear subcompo- 
nent is made by a starburst with metallicity ~ 2.5-6 times 
that of the bulge population and nuclear line-of-sight fb 
ranging from 0.13 to 0.91. Low /f,'s are needed for recent 
bursts and vice versa. Extreme /t,'s satisfy the constraints 
marginally whilst fb between about 0.33 & 0.91 give bet- 
ter fits. The burst mass may be lowered if its metallicity 
is raised, as discussed in § H. Age estimates range from 0.4 
Gyr to 1.7 Gyr, depending on fb, Zbuige & Znuc Tighter 
constraints would require better data, although the param- 
eters would somewhat be limited by degeneracy. 

7.1.3 Starbursts with enhanced [Fe] abundance (SB6) 

If the kinematic subcomponent is formed from gas accreted 
from a spiral into a giant elliptical, it may have a [Mg/Fe] 
abundance ratio lower than that of the underlying ellipti- 
cal, because giant ellipticals are thought to have an over- 
abundance in light-elements (or under-abundance in heavy 



7.2.1 



elements). This is further discussed in 

There are 2 ways of approach. The first is to keep the 
Fe abundances in the 2 galaxies the same, and raise the 
Mg in the elliptical. This will not work as too much depres- 
sion in <Mg> and Fe5335, and too little H/3 enhancement 
is predicted. The second way, which is adopted, is to keep 
the Mg abundances the same, but raise the Fe in the burst 
population. In this case Znuc is actually slightly higher than 
Zbuige- Because stellar populations in GISSEL96 have solar 
[Mg/Fe] abundances, we artificially enhance the burst Fe5335 
by decreasing the flux in model passbands, whilst ignoring 
the physical processes involved. Fig. lid show the results for 
model SB6 where Znuc ~ Zbuige = IZq, with the burst Fe 
enhanced by 0.08 dex over solar (model SB6). We find this 
model is marginally ruled out. Although it may be saved by 
increasing the Fe abundance by more than 0.08 de x, it is 
prob ably difficult to increase by much, as fig. 1 of (Davies 
199(:) shows that the typical abundance difference is ~ 0.1 



dex. 



7.1.4 Constraints on the burst duration (SB7, SB8) 

The burst duration is constrained by repeating model SB3 
with burst durations of 0.2 and 0.4 Gyr, respectively, in mod- 
els SB7 & SB8. The line-indices for the latter are plotted in 
Fig. |ll|. The observational constraints can be marginally met 
0.5 Gyr after the burst with /;, = 0.20, after 0.6 Gyr with 
fb = 0.33, or after 0.9 Gyr with fb = 0.91. The burst du- 
ration can be constrained weakly, and no solution exists for 
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Figure 8. The H/3, <Mg> and Fessss for model SBl where the nuclear subcomponent is formed in a 0.1 Gyr starburst with the same 
(solar) metallicity as the bulge. The ordinate is the line index and the abscissa is the galaxy age. The bulge is formed between & 1 Gyr 
and its line indices are plotted with thick solid lines. The nuclear component begins to form at 7 Gyr. Labelled along the right edges 
are the burst fraction /{, of the starburst. The observational constraints are shaded. In the central panel, the hatched region enclosing 
a dashed curve with square symbols indicates a nuclear <Mg> lowered by 0.053 ± 0.012 mag relative to the bulge. This same region 
(hatched) is mapped into the left and right panels for comparison purpose. In addition, the region corresponding to a nuclear H/3 enhanced 
by > 1.4A relative to the bulge is shaded in the left panel, whilst in the right panel, the region corresponding to lAFessssI < 0.35A is 
shaded. In order to have consistent solutions, the shaded and hatched regions must overlap in the left and right panels at the same age 
and /(,. 



durations longer than ~ 0.4 Gyr. Due to the short lifetimes 
of A- type stars there is less H/3 for a given /;,. As a conse- 
quence of increasing the burst duration, the observation has 
to be made earlier, and furthermore the range of permitted 
/b's becomes more restricted. 

We conclude that the burst duration can only be weakly 
constrained. Shorter bursts are favored because the range of 
/i, is less restricted. We also note that because the orbital 
period of the secondary galaxy is an order of magnitude 
lower than the maximum burst duration, it will probably be 
impossible to distinguish a succession of bursts from a single 
burst. 



7.2 Truncated star formation (TSF1-TSF3) 

Next we investigate whether the nuclear subcomponent is 
formed by continuous star formation that has been trun- 
cated recently. The feasibility of this model is further dis- 
cussed in § y. The stars in the nuclear subcomponent are 
formed between 1 & 7 Gyr with the same, solar metallicity 
as the bulge. The indices are plotted in Fig. [13 for model 
TSFl. The strength of star-formation along the line of sight 
is parameterized by /;,, redefined as the mass of the nuclear 
subcomponent divided by the total mass. The observational 
constraints cannot be met. Even with /t = 0.99, the desired 
AH/? is only attainable some 0.2 Gyr after tstop, whilst the 
AFe5335 of — 0.75A is ruled out at 98% confidence level. A 
high mass is required because only the recently formed stars 
are responsible for lowering the nuclear metallic line indices. 
As a consequence too many stars are made during prolonged 
star formation. This problem will be worse if the bulge has 
sub-solar metallicity so that the galaxy is older. 

This model cannot be saved by lowering or increasing 
the metallicity of the accreted stars either. If Znuc < Zbuige, 
the desired A<Mg> can be reached with a smaller ft, but 



there are not enough young stars to give the desired AH/3. 
In addition, AFessas will be too large. If Znuc > Zbuige, 
AFe5335 will be smaller, but an unrealistically high /;, is 
required to give the desired A<Mg>, because A<Mg> is 
pushed up too by the increased metallicity. This is confirmed 
when model TSFl is repeated with Z„uc = 2.5Z0, which 
shows that no solution exists for this model. 



7.2.1 Enhanced Fe abundance in the accreted stars 
(TSF2, TSF3) 

Giant ellipticals tend to have Mg indices enhanced relative 
to Fe (or Fe is depleted relative to Mg), whilst spiral galax- 
ies tend to have solar [Mg/Fe] abundance ratios. If the kine- 
matic subcomponent consist of stars accreted from a spiral 
into a giant elliptical, it is likely to have a [Mg/Fe] ratio lower 
than that of underlying elliptical. We find a difference in 
[Mg/Fe] abundance ratio alone can explain the observations. 
This is demonstrated when the above model is repeated with 
the Fe abundance of the nuclear subcomponent artificially 
enhanced by 0.04 & 0.08 dex respectively (models TSF2 & 
TSF3), whilst keeping the metallicity in both components 
solar. Model TSF2 is marginally ruled out by observations. 
The Fe5335 indices for model TSF3 are plotted in the right 
panel of Fig. h3, showing that an Fe enhancement of 0.08 
dex can meet the observational constraints, although rather 
high /b's are required. Solutions exist at At = 0.2 Gyr for 
0.67 < fb < 0.80, or at 0.4 Gyr for fb = 0.99. Therefore 
a starburst is unnecessary if the nuclear population has a 
Fe abundance enhanced by 0.08 dex, and are accreted by a 
process which truncated its star formation. A rather high 
fb ~ 0.50 is required because only the stars made recently 
are responsible for the enhanced H/?, thus many stars are 
made over a Hubble time. 
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Figure 9. The H/3, <Mg> and Fe5335 for models SB3— SB5 where the nuclear subcomponent is formed in a 0.1 Gyr burst with metallicity 
higher than that of the bulge. The top row corresponds to ^i,„;ge = I^q and Znuc = 2.5 Zq (model SB3), the middle row corresponds 



to ^buige = 0.4 Zq and Z„ 



1 Zq (model SB4), whilst the bottom row corresponds to Zi^^ig^ = 0.4 Zq and Znuc = 2.5 Zq (rnodel 



SB5). Labelled along the right edges are the burst fraction /;, of the starburst. The observational constraints are shaded; see Fig H for 
explanation. 



7.3 Continuous star formation up to the present 
(CSFl, CSF2) 

Is the nuclear subcomponent a primordial, undisturbed, 
spiral-like disk? In this scenario it is formed continuously 
from 1 Gyr onwards up to the present, and has the same 



metallicity and abundance as the bulge. For the purpose of 
this work, the bulge mass is arbitrary and set to IMq. Line- 
index predictions from the Big Bang onwards are presented 
in Fig. hsi for model CSFl where Z„uc = Ztuige = 1, with 
SFR parameterized in units of Mg/Gyr. This model cannot 
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Figure 10. The line indices for model SB7 where the nuclear component is formed in a 0.1 Gyr starburst at 7 Gyr, with the same (solar) 
metallicity as the bulge. However, the bulge Fe abuncance is 0.08 dex over the solar value (model SB7). Labelled along the right edges 
are the burst fraction /(, of the starburst. The observational constraints are shaded; please see Fig N for explanation. 
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Figure 11. The line indices for model SB8 where the nuclear component is formed in a starburst which started at 7 Gyr, and lasted 0.4 
Gyr. The burst has metallicity of 2.5 ZQwhilst the bulge has solar metallicity. 




Figure 12. Left 3 panels: The H/3, <Mg> and Fessss for model TSFl where the nuclear subcomponent is formed by constant star 
formation, truncated at 7 Gyr. The metallicity of the nuclear subcomponent is the same as that of the bulge, which is solar. Labelled 
along the right edges are the burst fraction /;, along the line-of-sight, defined as the mass of the nuclear component divided by the total 
mass. Right panel: The Fe5335 for model TSF3 in which the accreted stars have a Fe abundance of 0.08 dex over solar. The H/3 and 
<Mg> in this model remain the same as those in model TSFl. 
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Figure 13. The H/3, <Mg> and Fesaas for model CSFl where the nuclear subcomponent is formed by continuous star formation up 
to the present. Both the nuclear subcomponent and the bulge have solar inetallicity. The bulge (which has a mass of IMq) is formed 
between & 1 Gyr. The nuclear subcomponent is formed by constant star formation from 1 Gyr onwards. Labelled along the right edges 
in units of MQ/Gyr are the normalised star formation rates. 
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Figure 14. The H/3, <Mg> and Fessss for model RCSFl where the nuclear subcomponent is formed by star formation that begins 
recently and continues up to the present. The bulge, which has a mass of IMq, is formed between & 1 Gyr, whereas the nuclear 
subcomponent is formed by constant star formation from 7 Gyr onwards. In this model both the bulge and the nuclear subcomponent 
have solar metallicity. Labelled along the right edges in units of MQ/Gyr are the normalised star formation rates. 



explain simultaneously the observed weakness of <Mg> and 
the strength of Fesaas (right hand panel of Fig. |l3[). Likewise 
model CSF2, where Znuc ~ Ztuige = 0.4 Zq cannot explain 
these parameters. 

These models are only approximations because self- 
enrichment is not modelled in GISSEL96. If there is self- 
enrichment, the metal indices in FigJlSl will have steeper 
gradients, with H/3 remaining the same due to its weak de- 
pendence on metallicity. Whilst AFesaas is reduced, a very 
high SFR is required to attain the observed A<Mg> be- 
cause many young stars are needed to dilute the enhanced 
metal lines in the Ught of the old stars. As a rough esti- 
mate, model CSFl is repeated with Z„uc ~ 2.5 Zq. We 
find that, along the line of sight, the nuclear subcomponent 
has to be ~ 11.2 times more massive as the bulge. Such a 
massive disk will make the isophotes significantly disky, not 
observed in NGC 2865. Therefore we conclude the nuclear 



subcomponent is not a primordial embedded disk that has 
been forming stars up to the present. 



7.4 Recent, continuous star formation 

(RCSF1,RCSF2) 

Finally, was the nuclear subcomponent formed recently by 
continuously feeding gas of constant metallicity into the 
galaxy centre, for example, by th e return of tidal materi al 
after a major spiral-spiral merger (Hibbard & Mihos 1995)? 
In this model, the new stars starts to form at tstart and 
continues constantly up to the present; their metallicity can 
be the same or higher than that of bulge. The bulge mass 
is arbitrary and is set to 1 Mq. The model line-indices for 
model RCSFl where Zb^ige = Znuc = 1 Zq are plotted in 
Fig |l4[ The observed AH/3 & A<Mg> are attainable for 
SFR = 0.05-0.1 MQ/Gyr within 3 Gyr after tstart, but the 
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predicted AFessss = —0.75 A is ruled out at 98% confidence 
level. We find that model RCSF2 is similarly ruled out. The 
above models are also repeated with Zbuige = 0.4 Zq and 
Znuc = 0.4 or IZ0, and we find that they are also ruled out. 
Therefore we conclude the nuclear subcomponent was not 
formed by continuous inflow of cold gas with fixed metallic- 
ity on a Gyr timescale. 

7.5 Summary of Population Synthesis Models 

The results of this section are tabulated in Table 3. The 
observed line strength indices allow the following models: 

- Recent burst of star formation with enhanced 
metallicity. A recent burst of star formation works only if 
the resulting KDC has a metallicity higher than the galaxy 
bulge by a factor of about 2.5-6.3 (models SB3, SB4 & SB5). 
The burst duration can be as long as 0.4 Gyr (model SB8). 
These models do not fit the data if the KDC has the same 
metallicity as the bulge, even if [Mg/Fe] ratio of the KDC is 
lower. 

- Truncated star formation in the KDC, with a 
difference in abundance ratio. In this case the metallic- 
ity of the 2 populations are the same, but the Fe abundance 
in the KDC is enhanced relative to that of the bulge by 0.08 
dex (model TSF3) . The model does not fit the data without 
this enhancement. 

Continuous Star Formation and Recent Continuous 
Star Fo rmation models are completelv ruled out by the data. 



outer shells, their blue colours, and their chaotic distribution 
(Fort et al. 1986). The weak interaction model of decoupled 
core formation (Hau & Thomson 1994) can potentially ex- 
plain the high vja of NGC 2865 as due to the spinning up 
of the bulge by another galaxy. However, in this model the 
nuclear disk is primordial, in conflict with the observed line 
indices. It is also unlikely that a passing galaxy can transfer 
a large amount of orbital angular momentum over a period 
longer than 0.5 Gyr without being captured or substantially 
disrupted, as NGC 2865 has an extended massive dark halo 
(Schiminovich et al. 1995). Thus a purely interaction in- 
duced origin for the shells and KDC in NGC 2865 is ruled 
out. 

A gas rich merger and starburst origin for the KDC? 

The presence of a young population in the galaxy core is 
evidence for a late merger event. Not only is the population 
young, but it may be formed in a burst, a typical result of 
a spiral-spiral (SS) merger, or of the accretion of a gas-rich 
object by an elliptical. The distinct core can be formed by a 
starburst with metallicity 2.5-6 times that of the underlying 
galaxy, with 0.13 < ft < 0.91 and a short ( < 0.4 Gyr) 
duration. The best flt age, which ran ges from 0.4 to 1.7 Gyr , 
agrees with the age determined by Bica & AUoin (1987). 



It depends on /;, and metallicity, with more recent bursts 
favored for low /;, or high metallicity, and vice versa. 

The range of metallicity enhancement required is in 



good agreement with those expected by theory (Fritze- 



Alvensleben fc Gerhard 1994 [Milios fc Hernquist 1994] ) 



8 DISCUSSION 



A common merger origin for shells and KDC? 

The discovery of a kinematically distinct component at the 
heart of a shell galaxy provides additional evidence linking 
shell systems to mergers or accretion. Indeed, the distinct 
kinematics at the core of NGC 2865 require formation dy- 
namics decoupled from those of the rest of the system, and 
a merger or accretion provides a natural mechanism for such 
a process. 

This finding strengthens t he link betwee n shells and 
kinematically decoupled cores. Forbes (1992) found that 
galaxies selected for their KDC's have a high probability 
of having shells. The present paper has found that a galaxy 
selected for its shells has a kinematically distinct core. In 
a forthcoming paper (Balcells, Hau & Carter, in prepara- 
tion), we show that NGC 474, a shell lenticular, shows core 
kinematic decoupling as well. 



By considering a 1-zone closed boxed model, Fritze-v. Al- 
vensleben & Gerhard find that immediately ( £ 0.1 Gyr) 
after a collision of a spiral pair at 12 Gyr which consumed 
all their remaining gas, the global metallicity enhancement 
factors are about 1.5, 1.9 & 7.5 for Sa, Sb and Sc pairs re- 
spectively (Figs 12a & 12b of Fritze-v. Alvensleben & Ger- 
hard 1994). The metalhcity of NGC 2865 is roughly solar, 
similar to galaxies earlier than Sc. Taking both metallic- 
ity enhancement and absolute metallicity into account, if 
NGC 2865 is the remnant of the merger of 2 spirals of equal 
Hubble type, they are most likely to be of type between Sb 
& Sc (although see later discussion). 

The metallicity observed in the peculiar core of 
NGC 2865 is highe r, by large factors, than that of early 
type spiral bulges ( Balcells fc Peletier 1994 ) and that of 
HII regions of disks (Vila-Costas & Edmunds 1992), hence 



self-enrichment must have taken place, suggesting a burst 
duration longer than SN II feedback timescale. SPH sim- 
ulations of spiral-spiral (SS) merg ers which include star- 
formation (Mihos & Hernquist 1996) show that galaxy struc- 



The weak interaction model for shells (Thomson 1991) 
does not predict the existence of a kinematically distinct 
nucleus, as the shells are produced as a reaction to a graz- 
ing incidence passage; the presence of the rapidly rotating 
nuclear system in a shell galaxy is purely coincidental in 
this framework, unless the core is formed by mass transfer. 
Schiminovich et al. 1995 note the presence of a nearby gas- 
rich dwarf, but conclude that it is not related to NGC 2865. 
Furthermore because in this model the shells are in an one- 
arm density wave, they should follow the brightness of the 
underlying elliptical, in confiict with the existence of bright 



ture plays a dominant role in regulating starburst activi- 
ties. If the progenitors have massive central bulges, gas in- 
flows are strongest at the flnal stage of merging, and the 
resulting starburst is intense but short-lived (~0. 02-0. 03 
Gyr) , whereas if the progenitors are bulgeless (or have small 
bulges), the starburst is strongest after the first passage of 
the galaxies, with successive episodes of star formation af- 
ter each passage, and the final merger consumes the rest of 
the gas. Therefore successive generations of starbursts have 
higher and higher metallicity. Furthermore, the burst dura- 
tions in bulgeless encounters are much longer (~ 0.15 Gyr), 
thus allowing more time for self-enrichment in each star- 
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burst, 'therefore a late type spiral progenitor with a small 



bulge is favored over an early type spiral with a large bulge. 
This will also be consistent with the shell colour indices 
which Fort et al. (1986) found to resemble those of an Sb or 
Sc galaxy. 

Although the above discussions are based on simula- 
tions of SS mergers, because gas inflow is driven by the tidal 
torque of the galaxy the gas resides in, the general picture 
above will probably be true for a elliptical-spiral (Es or ES) 
merger too. 



The identity and mass ratio of the progenitors 

Does our data allow us to clarify the discrepant views on the 
type of merger which formed NGC 2865 (§ 1) ? If the s hells 
contribute ~ 11 — 22% of the total light of the galaxy (Fort 



et al. IS 186), then taking into account that only a small frac- 
tion of a merged galaxy ends up in visible shells, the fraction 
of the total system mass intially in dynamically cold compo- 
nents cannot be much lower than 0.5. Hence, if one galaxy 
was an elliptical, the galaxy mass ratio must have been 
close to unity. This argues against the minor merger scenario 
proposed by Bettoni (1992). Alternatively, both progenitors 
could have been disk galaxies. In the disk-disk case, we can 
constrain the mass ratio to below a value around 3:1 or 5:1, 
as high er mass ratios do not destroy the disk of the larger 
galaxy (Barnes 1996). The large amount of rotational sup- 



port {v/g ~ 1.2), one of the highest known for an elliptical, 
could be a relic of the disk rotation of the larger precursor in 
an intermediate mass-ratio merger. It could also result from 
a good alignment of the two spins in an 1:1 SS merger. The 
good alignment between stellar and gas kinematics could be 
a relic of disk rotation in an intermediate mass merger of two 
spirals, or the result of an organizing effect of the merger in 
an equal mass SS merger. 



Discrepancy between age of shells and the KDC- against a 
phase-wrapping origin of shells? 

The ages estimated for the NGC 2865 starburst are higher 
than ty pical ages for shell system s. Using phase wrapping 
theory, Hernquist & Quinn (1987) derive an age of ~ 5 x 
" WhT,l, yr for NGC 1344 and 



noo 



10 ftroo yr for NGC 3923, 

~ 2 X 10® /i]"qq yr for NGC 3051. For a given mass model, 

the tim^ since tne merger scales as 



Tn 



(n-f0.5)P(d„) 



(6) 



(formula [5.3] of Hernquist & Quinn (1987)), where n is the 
number of shells and P{d„) is the time between turning 
points (half the orbital period) at the location of shell n 
(the innermost shell) . Large galaxy sizes and a high number 
of shells translate into a long merger age. Notwithstanding 
that the chaotic nature of the shells in NGC 2865 precludes 
a rigorous application of phase-wrapping theory, we apply 
this formula directly, seeking an order-of-magnitude result. 
We further approximate 



P{dn) « 2d„/a, 



(7) 



which should be acceptable for our purposes if d„ « re. We 
take n = 9, i.e. the 7 shells identified b y Fort et al. (1986) 
plus the two outer loops seen in Fig. lb of Schiminovich et al 



2-1 ''■100 kpc 



(1995|); dn = 20" = 2.5 fe7nn kpc, r^ = 16.6" = 

( p0rgensen, Franx, fc Kjaergaard 1992 ), and a = 200 kms~^ 

(this paper). We obtain 



in 



2.4 X 10*^ h^olyr 



(8) 



Hence, phase- wrapping yields tmerger an order of magnitude 
below the burst age derived from population synthesis. Un- 
certainties in the true number of shells, or in the number 
of periods completed by stars in the outermost shell, can at 
most introduce a relative error of a few tenths in tmerger, as 
they contribute additive terms of order 1 to n in formula pi 
Approximation (frl) for P{d„), when applied to NGC 3923, 
gives t 



merger 



in accord with the full model calculation of 



Hernquist & Quinn (1987) to within 25%. Thus, it is un- 



likely that the approximations in the calculation can account 
for the discrepancy between the phase-wrapping merger age 
and the age derived for the stellar population. In theory, 
the age of the KDC can be shortened by raising its metal- 
licity to higher than 6.3 times that of the the bulge. More 
likely, this discrepancy indicates that the shell system did 
not form by phase- wrapping of a dynamically-cold, accreted 
satellite, and that the phase wrapping theory is not appli- 
cable to NGC 2865, either because of dynamical friction or 
because a major, rather than a minor, merger took place. 
Alternative mechanisms include spatial wrapping or phase 
wrapping in a SS merger remnant. These are supported by 
Fort et al. (1986), who remarked on the similarity between 
NGC 2865 and parabolic-encounter models (e.g. model 8 in 
Hernquist & Quinn 1987), where shells are formed by spa- 
tial wrapping over several orbits, whose formation timescale 
could be as long as 4 Gyr. 

Truncated star formation? 

Another way out of the age discrepancy problem is to form 
the KDC by accreting stars from a small star-forming disk 
galaxy into a giant elliptical, during which the star formation 
in the former is truncated. The two progenitor galaxies have 
the same Mg abundance, but the accreted population has a 
relatively higher Fe abundance. A high /;, of 0.5 ^ ft ~ 0.99 
is required. The age estimates of 0.1 - 0.4 Gyr is in bet- 
ter agreement with the age of the phase-wrapped shells. 
The physical processes involved in this model are contrived. 
An accreted stellar system can create a central co ncentra- 
tion in the rem nant via dynamical friction only (Balcells 



& Quinn 1990), and a small spiral is likely to tidally dis- 



rupt before reaching the center of the elliptical. However, 
it is conceivable that a relatively gas-poor and small spiral 
sinks undisrupted thanks to the self gravity of its dark halo, 
forming the kinematically distinct core which dominates the 
nuclear light, whilst tidally-ejecting parts of its disk to be- 
come phase-wrapped shells. We note that if only the bulge 
remains bound, however, there will be a problem matching 
the age and metallicity observed in the center, as bulges of 
small spirals are not metal rich (Balcells & Peletier 1994). 



It is also unclear how to truncate the star formation in 
a star-forming galaxy during a merger. Below, we outline 
briefly the possibilities — quantitative discussion is beyond 
the scope of this work. It may be possible that during the 
merger process, young stars are formed by tidal shocking 
of gas along transient spiral arms, but they are not con- 
flned to the nucleus, where the old spiral disk dominates the 
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light. Another possibihty is the evaporation of the spiral HI 
by the hot ISM of the elliptical, which may effectively and 
cleanly truncate the star formation, but may only work for 
a very small accreted galaxy. We note that SPH models o f 
elliptical-spiral (Es) mergers (e.g. Hernquist & Weil 1992) 



Mg^ 



do not model the hot ISM in the elliptical. Other possibil- 
ities are ram-pressure stripping of HI by the hot elliptical 
ISM, or the stripping of the remaining HI by the winds of 
the first supernovae. Ram-pressure stripping could possibly 
ex plain the displacement of HI outside the shells reported 
by Schiminovich et al. (1995 ). 

Finally, given that both a starburst and a truncated 
star-formation scenario explain the line indices, it is possible 
that both processes played a role in NGC 2865. Indeed, a 
natural star formation history for a gas-rich galaxy may be 
that of forming stars throughout its life, and at some point 
having its remaining gas consumed in a starburst after a 
collision. 



Can the age difference between the shells and the KDC 
reflect the merging time? 

A further possibility is that the shells and the decoupled 
core are not formed at the same time. This is quite plausi- 
ble, as phase- wrapped shells are formed at the final stage of a 
merger, whereas starburst activities are regulated by the in- 
ternal structure of the spiral, and maybe strongest after the 
first pa ssage and well before the final merger for a bulgeless 
spiral (Mihos & Hernquist 1996). Therefore if the shells are 



indeed phase-wrapped, the age difference between them and 
the KDC may refiect the merging time, and further favor a 
late-type spiral progenitor. 



A comment on Mg enrichment and the adopted IMF 

Throughout this paper, a Scalo IMF has been adopted. A 
proper investigation on the choice of IMF is beyond the 
scope of this work, and is further hampered by the fixed 
abundance ratios in the models, and the lack of knowledge 
on self-enrichment of individual elements. Starbursts are sus- 
pected to produce more massive stars, which are thought to 
feedback more Mg than Fe into the ISM through type II 
supernovae. We find in § o that the <Mg> of a burst popu- 
lation is weaker than the <Fe> before 1.7 Gyr, due to a tem- 
perature effect in the atmospheres of stars at the MS turnoff. 
This helps to hide any enhanced Mg abundance during this 
time. Therefore, unless the Mg is much more enriched than 
the Fe during the starburst, our results are probably insen- 
sitive to the Mg enrichment. This work would benefit from a 
better understanding of the role SN feedback, gas dynamics 
and star-formation play in gas-rich mergers. 



The 
Mg- 



'hulge+burst" assumption, and the departure from the 
c( relationship 



The very "old" nature of the underlying population is an 
assumption of the synthesis models. Outside the young core 
(r > 10") H/3 « 2A (Fig. §). A roughly 10% pollution by 
stars of order 1 Gyr old onto an u nderlying 12 Gvr popula - 
tion results in such an H/3 value (De Jong fc Davies 1997). 
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Figure 15. The effect of light contamination by young stars on 
an old bulge population. The bulge is formed between & 1 Gyr 
with solar metallicity. Its Mg2 is indicated by the solid line. Added 
to this population are young stars of the same metallicity formed 
between 12.0 and 12.1 Gyr. The dashed line shows the Mg2 for 
the combined light, where the young stars contribute 9% of the 
total mass. The dotted line shows the same for a 17% contribution 
by young stars. 



This will also be consistent with the blue shells which con- 
tribute ~ 11 — 22% of the total luminosity. In § Wwe find that 
the Mg2 over the galaxy is about 0.08 mag below that ex- 
pected from the Mg — a relationship. This has also been 



noted independently by j0rgensen (1997). Is this scatter 



from the fundamental plane due to light contamination by a 
young population? To investigate this, we add young stars to 
an old solar-metallicity population, and the model Mg2 are 
plotted in Fig. [13. We find that a roughly 10% contamina- 
tion by 1 Gyr old stars can indeed result in the magnitude of 
scatter observed. Therefore the young population is not only 
confined to the nucleus, but a fraction is also mixed into the 
underlying elliptical population. This makes the galaxy ap- 
pear younger or has a lower metallicity. In light of this, the 
models in table 3 with a solar metallicity bulge are probably 
more relevant than those with sub-solar metallicity. 

NGC 2865 shows that gas-rich, shell-forming mergers 
can be a source of scatter from the Mg — a relationship. 
The "bulge-|-burst" assumption is an approximation which 
helps to keep the models simple. But because the light con- 
tamination by the young stars on the bulge is small, the 
conclusions from population synthesis are unlikely to be sig- 
nificantly affected. 



KDC formation issues 

Our discovery of a KDC in a shell elliptical galaxy, with as- 
sociated young stellar population, fits well in the framework 
of gas rich spiral merger sim ulations (Hernquist & Barnes 



1991| , [Hernquist fc Weil 19921 ). The fact that the light profile 
of NGC 2865 is well fit by a r 3 law, with a kinematically 
distinct, compact and young component embedded in a dy- 
namically hot system, is consistent with the segregation of 
stars and gas seen in these simulations. 

We demonstrated that despite the age and metallic- 
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ity degeneracy of the underlying elliptical, physical param- 
eters such as the age and metallicity of the kinematic sub- 
component in a recent merger can be constrained. In addi- 
tion we discovered that, between about 0.4-1.8 Gyr after a 
metallicity-enhanced starburst, the <Fe> of stars at the MS 
turnoff is strong enough to compensate for dilution effects, 
whilst <Mg> is weaker. As a result, the abundance ratio 
of the total light appears to be modified. This is purely a 
temperature effect in the atmospheres of stars at the MS 
turnoff. As a consequence of the weaker Mg index in the 
burst population, Mg enrichment by a starburst may be par- 
tially hidden during approximately 1.8 Gyr. 

In a sample of galaxies chosen subjectively by their vi- 
sual appearance to represent the time sequence of SS merg- 
ers, NGC 2865 was ranked as one of the oldest remnant 



identifiable by its morphological features ( Keel fc Wu 1995 ) 



That ranking was based on the basis of the shells and de- 
viations from elliptical shape. The dating of such a galaxy 
by morphology is difficult because tidal tail remnants dis- 
appear after a few crossing times. It is also unlikely that 
visual appearance alone can provide enough evidence for 
a SS merge r origin for a r emnant like NG C 2865, unlik e 
NGC 7252 dSchweizer 198^ or NGC 3656 ( [Balcells 199?! ). 
This work demonstrates that an older merger remnant may 
still be dated by its stellar populations, even when the tidal 
tail has virtually disappeared. 

NGC 2865 has the only KDC with depressed Mg 
indices- most, if not all, decoupled cores have enhanced Mg , 



in c omm on with most giant ellipticals (Bender & Surma 
1992 |). p ointing to the different [Mg/Fe] abundance ratios 
between spirals and giant ell ipticals, and t he fact that KDCs 
tend to reside in the latter. Bender (1996) favored the view 
that KDCs are metallicity-enriched remnants of early hi- 
erarchical mergers rather than remnants of present-day S- 
S mergers. We find that NGC 2865 could be the renmant 
of a recent S-S or S-E merger. What would the metallicity 
of the decoupled core be in the next few Gyrs? If it was 
formed by a metallicity enhanced starburst, then its metal 
line-indices would eventually be higher than those of the 
underlying elliptical. If the IMF is top-heavy and the star- 
formation timescale is short, there would be an additional 
Mg en hancem ent observable after ~ 1.9 Gyr ( Matteucci & 
Brocatd 199C). Thus a recent starburst origin will fit the 



overall scheme of KDC formation. NGC 2865 is an example 
of a young KDC fo rmed in a recent merge r. Unlike for in- 
stance NGC 7626 ( [Balcells fc Carter 1993| ) NGC 2865 has 
symmetric, regular velocity field. It is unclear what propor- 
tion of young KDCs do have a regular velocity field, and 
therefore how quickly the velocity field relaxes. 
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APPENDIX A: RECOVERABILITY OF 
KINEMATIC PARAMETERS IN THE 
PRESENCE OF NOISE 

The sensitivity of the kinematic parameters to noise is inves- 
tigated by simulating the effect of both Poisson shot noise 
and residuals from imperfect sky subtraction. The artifi- 
cial data is generated by convolving a stellar template with 
Gaussian and non-Gaussian LOSVDs of cz = 2630 kms~^ 
and a = 190 kms^^, scaled to different continuum levels, 
and added Poisson noise to make spectra with S/N — 40, 
55 and 70, typical of our data. The spectra at the nucleus of 
NGC 2865 have S/N ~ 70, whilst those at intermediate dis- 
tances (~ 2 -3 ) have S/N ~ 40. The effect of imperfect sky 
subtraction is mimicked by adding rows of real 'sky' spectra 
extracted sufficiently away from the galaxy. For a Gaussian 
LOSVD, up to 12 rows are added to mimic the worst case 
of sky contamination where heavy summation is required 
to attain the desired S/N ( > 40) typical of the outermost 
points. Recovery of non-Gaussian LOSVDs is tested by set- 
ting one Hermite moment to 0.05 whilst keeping the others 
zero. The kinematic parameters are recovered 3 times us- 
ing kinematics, each with different Poisson noise seeds and 
sections of blank sky. 



Perfect Gaussian h,=0.05 h,=0.05 h..=0.05 hj.=0.05 



The results of the investigation are plotted in Fig. Al 
The mean velocity v can be determined to within ^ 3 kms~ 
at the galaxy nucleus, and ~ 20kms~^ at the outermost 
points, ft is insensitive to asymmetric deviations {h^ & h(,), 
with only a small positive offset when h^ is introduced to 
the LOSVD, and no significant offset when hs is introduced. 

If the LOSVD is Gaussian, the velocity dispersion a 
is at best determined to within ~ 3kms~^ at the nucleus. 
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Figure Al. The effect of noise to the recovery of the parame- 
ters. The composite stellar template at the nucleus of NGC 2865 
is convolved with different model LOSVDs to make mock galaxy 
spectra, and added to them Poisson noise and blank sky spectra. 
The kinematic parameters are recovered 3 times, each with a dif- 
ferent noise seed and sections of blank sky. The leftmost panel 
corresponds to results for a Gaussian LOSVD, and the next 4 
panels correspond to results with a signal of 0.05 in each of h^, 
/i4, /15 &; hg respectively. The model values are represented by 
the dotted lines, and the recovered parameter the solid dots. In 
each panel, results are shown for S/N of 40, 55 and 70, labelled 
along the bottom with the number of rows of sky summed. For 
a Gaussian LOSVD, 5 and 12 rows of sky spectra are added to 
the artificial spectrum to mimmick the points farthest from the 
galaxy nucleus. 



and to -^ 5 kms~^ at intermediate distances. Non-Gaussian 
LOSVDs increase the scatter further, with the introduction 
of 0.05 in each Hermite moment increasing the scatter in 
a by about 2kms~^ at the nucleus, and about 5kms~^ at 
intermediate distances. In addition, an /14 of 0.05 introduces 
a systematic displacement of ~ +4 kms~^ in a. 

In general, the Hermite moments h^, h^, h^ and /ig can 
be recovered with reasonable success, although with decreas- 
ing S/N , the signal cannot be fully recovered. Furthermore, 
signals in higher-order Hermite moments tend to leak into 
the moment 2 orders below. This effect is independent to 
S/N, and is probably a consequence of the finite resolution 
of the data. A Hermite moment h^ of 0.05 can be recovered 
to within 0.01 at the galaxy nucleus, and to within 0.02 at in- 
termediate distances. An hz, of 0.05 introduces a systematic 
displacement in /13 of about -1-0.01. Apart from a small pos- 
itive displacement in velocity, no systematic displacement is 
introduced in the other parameters when /13 is set to 0.05. 
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An hi of 0.05 can be recovered to within 0.01 at the nucleus, 
and within about 0.02 at intermediate distances. Non-zero 
values in other Hermite moments tend to increase the scat- 
ter in hi slightly. In addition, a systematic displacement of 
~ -1-0.02 in hi is introduced when he, is 0.05. 

The higher moments /15 and h% can be recovered reason- 
ably well, with a scatter of 0.01 at the nucleus and 0.03 for 
the intermediate points. The h^ cannot be recovered fully, 
and is systematically lowered by about 0.01. 

There appears to be little problems associated with 
skyline-residuals from imperfect sky-subtraction. Even when 
12 rows of spectra are combined to attain the desirable S/N , 
the scatter in the recovered parameters are increased only 
slightly compared with those derived from high-quality spec- 
tra, and only with v most affected. 

The results from these exercises show that Hermite mo- 
ments up to /i6 can be recovered reasonable well for the 
points inside ~ 3 . The signals in the lower Hermite mo- 
ments may be affected if there is significant signal in the 
Hermite moment 2 orders up. The parameter least sensitive 
to noise is the mean velocity, whilst the velocity dispersion 
G is affected by lowering the S/N or by increasing the power 
in the Hermite moments. 



APPENDIX B: DEPENDENCE OF THE 
KINEMATIC PARAMETERS ON THE CHOICE 
OF TEMPLATES 

Apart from noise, the kinematic parameters are sensitive to 
template mismatch, which could be caused by a difference 
in stellar type between the galaxy population and the tem- 
plate, or by a difference in metallicity abundance. The for- 
mer can be somewhat minimised by utilisation of an optimal 
template, taken as the best linear combination of the stellar 
templates (spanning a range of stellar type) that minimises 
the x^5 after an initial estimation of the LOSVD. In theory, 
the latter could be minimised too if the abundances of the 
template stars span a range that brackets that of the galaxy. 
The Lick indices of the template stars span the range 0.015- 
0.096 mag in Mgi, 0.090-0.220 mag in Mg2, 1.52-3.43 A in 
Mg 6, 3.07-4.30 A in Fe5270 and 2.74-3.92 in Fessss. Thus 
the Mgi , Mg2 and Fessss are well matched, and the Mg b and 
Fes270 marginally matched. Thus template mismatch due to 
metal abundance difference is unlikely to be a problem. We 
note if absorption lines of the stars are much weaker than 
that of the galaxy, then the hi may be artificially enhanced 
as the model spectrum is made more "pointy" to match the 
galaxy spectrum. 

For the purpose of quantifying the sensitivity of the 
Gauss-Hermite parameters to template mismatch, they are 
derived separately with each of the 3 available stellar tem- 
plates. The results are plotted in Fig. A.2. The parameters's 



Three different template' 



zero points vary with the choice of template but their ra- 
dial profiles do not. The extreme differences are 6kms~^ 
for both « & cr, and 0.026, 0.026, 0.020 and 0.018 for hz, 
hi, hs & he respectively. These should be regarded as the 
maximum zero point offsets, because the employment of an 
optimal template ensures that template mismatch is mini- 
mized. We conclude that the radial profiles of the Gauss- 
Hermite moments, apart from hi, cannot be explained by 
radial template mismatch even at the most extreme. The 
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Figure A2. The Gauss-Hermite parameters he, /is, hi, /is, a 
and V recovered separately using the 3 available stellar templates, 
showing that their zero-points depend on the choice of template, 
but not their radial profile. 



lowered point in /14 at r = 0.3 is reproduceable by extreme 
template mismatch. 

In the final analysis (see fig. H of main text), the zero- 
point of ~ -1-0.02 in /is is explainable by mismatch, and 
likewise /is and he- The absolute zero-point of the /14 is re- 
lated to the radial anisotropy of the orbital population. The 
constant /14 offset of ~ -1-0.045 along the major axis is larger 
than the extreme difference of 0.026 introduceable by the 3 
templates. Even taking into account the result from the pre- 
vious section where /ig ~ 0.02 may introduce an additional 
offset of ~ 0.01 into the /14, its offset is still marginally sig- 
nificant. 
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Table Al. Measured parameters along the major axis of NGC 2865. 
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2614.0 


1.9 


175.6 


1.8 


0.028 


0.013 


0.034 


0.012 


0.016 


0.014 


0.016 


0.014 


1.3 


2595.7 


2.6 


180.2 


2.1 


0.041 


0.016 


0.028 


0.016 


0.020 


0.017 


0.018 


0.017 


2.0 


2588.1 


3.7 


187.1 


3.2 


0.058 


0.016 


0.025 


0.016 


0.012 


0.017 


0.014 


0.017 


3.2 


2599.3 


3.8 


192.3 


3.5 


0.059 


0.016 


0.022 


0.016 


0.018 


0.017 


0.046 


0.017 


5.9 


2588.0 


3.2 


197.6 


3.3 


0.024 


0.015 


0.030 


0.015 


0.012 


0.016 


- 


- 


12.4 


2564.8 


3.6 


183.4 


3.0 


0.038 


0.016 


0.011 


0.015 


- 


- 


- 


- 



Table A2. Measured parameters along the minor axis of NGC 2865. 



K") 


V 


erry 


fj 


errtj 


hi 


errh3 


hi 


err^^ 


ht, 


err^j. 


he 


errhg 


-4.8 


2632.4 


3.7 


197.0 


3.6 


0.017 


0.017 


-0.007 


0.017 


_ 


_ 


_ 


_ 


-2.5 


2640.5 


3.1 


205.7 


3.1 


0.020 


0.014 


0.024 


0.014 


0.011 


0.014 


- 


- 


-1.4 


2632.7 


2.6 


184.8 


2.4 


0.010 


0.017 


0.036 


0.016 


0.037 


0.018 


0.039 


0.014 


-0.6 


2629.1 


2.4 


185.5 


1.5 


-0.008 


0.013 


0.016 


0.013 


0.029 


0.013 


0.014 


0.013 


0.2 


2630.9 


1.7 


192.8 


1.6 


0.002 


0.012 


0.016 


0.012 


0.019 


0.012 


0.023 


0.013 


1.0 


2638.7 


2.4 


185.5 


2.2 


-0.033 


0.016 


0.026 


0.015 


0.023 


0.017 


0.031 


0.017 


1.7 


2629.9 


3.5 


190.6 


3.2 


0.016 


0.016 


0.035 


0.015 


0.018 


0.017 


0.050 


0.017 


2.9 


2629.9 


3.6 


188.3 


3.5 


0.018 


0.017 


0.046 


0.016 


0.067 


0.018 


- 


- 


6.0 


2624.1 


4.0 


202.2 


3.5 


0.008 


0.016 


0.052 


0.016 


- 


- 


- 


- 


12.5 


2615.2 


5.1 


175.3 


4.4 


0.012 


0.022 


- 


- 


- 


- 


- 


- 
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Table A3. Measured line indices and their errors along the major axis of NGC 2865. 



r(") 


Hp 


err 


Mgi 


err 


Mga 


err 


Mg b 


err 


Fe5270 


err 


^65335 


err 


-14.92 


_ 


_ 


0.070 


0.016 


0.180 


0.020 


4.07 


0.72 


_ 


_ 


_ 


_ 


-14.15 


2.57 


0.37 


0.070 


0.016 


0.184 


0.020 


2.46 


0.73 


2.72 


0.43 


3.41 


0.41 


-13.38 


- 


- 


0.071 


0.015 


0.208 


0.019 


3.33 


0.69 


- 


- 


- 


- 


-12.61 


- 


- 


0.067 


0.014 


0.198 


0.018 


2.35 


0.68 


- 


- 


- 


- 


-11.84 


2.41 


0.33 


0.074 


0.014 


0.188 


0.018 


4.52 


0.61 


2.91 


0.38 


2.39 


0.38 


-11.07 


- 


- 


0.067 


0.013 


0.194 


0.017 


3.42 


0.61 


- 


- 


- 


- 


-10.30 


3.46 


0.51 


0.068 


0.013 


0.190 


0.016 


3.88 


0.57 


2.23 


0.61 


2.84 


0.60 


-9.53 


1.74 


0.50 


0.076 


0.012 


0.194 


0.015 


3.63 


0.55 


2.67 


0.57 


2.68 


0.57 


-8.76 


2.82 


0.47 


0.074 


0.011 


0.188 


0.014 


3.60 


0.51 


2.41 


0.54 


2.53 


0.53 


-7.99 


1.77 


0.45 


0.066 


0.011 


0.193 


0.013 


4.24 


0.48 


3.06 


0.51 


2.89 


0.50 


-7.22 


2.07 


0.42 


0.078 


0.010 


0.199 


0.013 


3.89 


0.45 


2.60 


0.48 


3.16 


0.47 


-6.45 


2.23 


0.39 


0.081 


0.009 


0.194 


0.012 


3.85 


0.42 


2.27 


0.45 


2.79 


0.44 


-5.68 


2.45 


0.37 


0.077 


0.009 


0.204 


0.011 


3.58 


0.39 


2.87 


0.41 


2.54 


0.41 


-4.91 


2.50 


0.33 


0.082 


0.008 


0.199 


0.010 


3.56 


0.36 


2.52 


0.38 


2.60 


0.37 


-4.14 


2.74 


0.30 


0.085 


0.007 


0.211 


0.009 


3.93 


0.32 


2.86 


0.34 


2.67 


0.34 


-3.37 


2.76 


0.27 


0.084 


0.006 


0.208 


0.008 


3.82 


0.29 


3.26 


0.30 


2.76 


0.30 


-2.60 


2.72 


0.22 


0.088 


0.005 


0.198 


0.006 


3.77 


0.23 


2.65 


0.25 


2.90 


0.24 


-1.83 


2.92 


0.18 


0.094 


0.004 


0.199 


0.005 


3.38 


0.18 


2.71 


0.19 


2.85 


0.18 


-1.06 


3.50 


0.15 


0.082 


0.003 


0.189 


0.004 


3.28 


0.13 


2.83 


0.13 


2.68 


0.13 


-0.29 


3.44 


0.13 


0.058 


0.002 


0.170 


0.003 


3.36 


0.11 


2.74 


0.11 


2.85 


0.11 


0.48 


3.45 


0.13 


0.037 


0.002 


0.153 


0.003 


3.31 


0.11 


2.82 


0.11 


2.76 


0.11 


1.25 


3.33 


0.14 


0.058 


0.003 


0.172 


0.004 


3.41 


0.14 


2.87 


0.14 


2.71 


0.14 


2.02 


3.29 


0.18 


0.076 


0.004 


0.192 


0.005 


3.51 


0.19 


2.87 


0.20 


2.75 


0.20 


2.79 


2.81 


0.22 


0.079 


0.005 


0.194 


0.007 


3.65 


0.24 


2.95 


0.25 


2.72 


0.25 


3.56 


2.73 


0.26 


0.075 


0.006 


0.197 


0.008 


3.93 


0.29 


2.61 


0.31 


2.73 


0.30 


4.33 


2.94 


0.30 


0.069 


0.007 


0.193 


0.009 


3.79 


0.33 


2.81 


0.35 


2.59 


0.34 


5.10 


3.20 


0.33 


0.063 


0.008 


0.202 


0.010 


4.15 


0.36 


2.84 


0.38 


2.82 


0.38 


5.87 


2.18 


0.36 


0.068 


0.009 


0.199 


0.011 


3.99 


0.39 


2.72 


0.42 


2.52 


0.42 


6.64 


2.33 


0.39 


0.069 


0.009 


0.173 


0.012 


3.71 


0.43 


2.63 


0.45 


3.15 


0.44 


7.41 


2.18 


0.42 


0.083 


0.010 


0.192 


0.013 


3.65 


0.46 


2.18 


0.49 


2.77 


0.48 


8.18 


2.41 


0.45 


0.073 


0.011 


0.196 


0.014 


4.27 


0.48 


2.18 


0.53 


2.66 


0.51 


8.95 


2.53 


0.47 


0.067 


0.012 


0.187 


0.014 


4.09 


0.51 


2.99 


0.55 


2.01 


0.56 


9.72 


1.89 


0.51 


0.077 


0.012 


0.184 


0.015 


3.68 


0.55 


2.40 


0.59 


2.76 


0.57 


10.49 


1.65 


0.54 


0.060 


0.013 


0.194 


0.016 


3.34 


0.59 


3.03 


0.60 


3.73 


0.59 


11.26 


- 


- 


0.051 


0.013 


0.177 


0.017 


4.23 


0.59 


- 


- 


- 


- 


12.03 


1.77 


0.34 


0.062 


0.014 


0.187 


0.018 


3.45 


0.63 


2.28 


0.38 


2.10 


0.38 


12.80 


- 


- 


0.058 


0.015 


0.161 


0.018 


3.63 


0.64 


- 


- 


- 


- 


13.57 


- 


- 


0.061 


0.015 


0.170 


0.019 


3.86 


0.67 


- 


- 


- 


- 


14.34 


2.24 


0.38 


0.047 


0.016 


0.178 


0.020 


3.77 


0.69 


2.34 


0.43 


2.45 


0.42 


15.11 


- 


- 


0.051 


0.017 


0.179 


0.021 


3.48 


0.73 


- 


- 


- 


- 
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Table A4. Measured line indices and their errors along the minor axis of NGC 2865. 



r(') 


^/3 


err 


Mgi 


err 


Mg2 


err 


Mg b 


err 


-^65270 


err 


^^65335 


err 


-5.20 


1.86 


0.36 


0.073 


0.009 


0.172 


0.011 


3.36 


0.39 


2.21 


0.42 


2.26 


0.41 


-4.43 


2.59 


0.33 


0.078 


0.008 


0.189 


0.010 


3.78 


0.36 


2.92 


0.38 


2.11 


0.38 


-3.66 


2.46 


0.30 


0.075 


0.007 


0.186 


0.009 


3.55 


0.32 


3.20 


0.33 


2.48 


0.34 


-2.89 


2.76 


0.25 


0.081 


0.006 


0.196 


0.007 


3.64 


0.27 


3.11 


0.29 


2.98 


0.28 


-2.12 


2.93 


0.21 


0.096 


0.005 


0.211 


0.006 


3.66 


0.21 


3.21 


0.22 


2.74 


0.22 


-1.35 


3.18 


0.16 


0.099 


0.003 


0.212 


0.004 


3.47 


0.15 


3.23 


0.15 


2.86 


0.15 


-0.58 


3.45 


0.13 


0.070 


0.002 


0.176 


0.003 


3.33 


0.11 


2.68 


0.11 


2.80 


0.11 


0.19 


3.49 


0.13 


0.043 


0.002 


0.162 


0.003 


3.43 


0.10 


2.82 


0.11 


2.91 


0.11 


0.96 


3.36 


0.14 


0.058 


0.003 


0.173 


0.004 


3.31 


0.13 


2.91 


0.14 


2.86 


0.14 


1.73 


3.20 


0.17 


0.077 


0.004 


0.191 


0.005 


3.52 


0.19 


2.81 


0.20 


2.69 


0.20 


2.50 


2.84 


0.22 


0.081 


0.005 


0.192 


0.007 


3.50 


0.25 


2.90 


0.26 


2.79 


0.26 


3.27 


2.53 


0.27 


0.066 


0.007 


0.193 


0.008 


3.89 


0.30 


2.87 


0.32 


2.38 


0.32 


4.04 


2.74 


0.31 


0.059 


0.008 


0.193 


0.010 


3.87 


0.35 


2.77 


0.37 


2.57 


0.37 


4.81 


2.09 


0.36 


0.076 


0.009 


0.182 


0.011 


3.96 


0.39 


2.74 


0.42 


2.83 


0.42 


5.58 


2.41 


0.40 


0.065 


0.010 


0.186 


0.012 


4.16 


0.44 


3.32 


0.47 


3.53 


0.46 


6.35 


1.93 


0.45 


0.073 


0.011 


0.196 


0.014 


3.27 


0.51 


2.96 


0.53 


3.29 


0.52 


7.12 


- 


- 


0.059 


0.012 


0.187 


0.015 


- 


- 


- 


- 


- 


- 


7.89 


2.29 


0.31 


0.089 


0.013 


0.181 


0.016 


4.24 


0.33 


2.65 


0.36 


2.10 


0.36 


8.66 


- 


- 


0.080 


0.014 


0.197 


0.017 


- 


- 


- 


- 


- 


- 


9.43 


- 


- 


0.072 


0.015 


0.140 


0.018 


- 


- 


- 


- 


- 


- 


10.20 


2.33 


0.37 


0.036 


0.016 


0.160 


0.020 


3.74 


0.40 


2.04 


0.44 


1.89 


0.44 


10.97 


- 


- 


0.057 


0.016 


0.233 


0.021 


- 


- 


- 


- 


- 


- 



